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O-i. ABSTRACT 

String and M-theory compactifications generically have compact mod- 
uli which can potentially act as the QCD axion. However, as demonstrated 
here, such a compact modulus can not play the role of a QCD axion and 
5^ i solve the strong CP problem if gravitational waves interpreted as aris- 

ing from inflation with Hubble constant Hi^f > 10^^ GeV are observed 
by the PLANCK polarimetry experiment. In this case axion fluctuations 
generated during inflation would leave a measurable isocurvature and/or 
non-Gaussian imprint in the spectrum of primordial temperature fluctua- 
tions. This conclusion is independent of any assumptions about the initial 
axion misalignment angle, how much of the dark matter is relic axions, 
or possible entropy release by a late decaying particle such as the saxion; 
it relies only on the mild assumption that the Peccei-Quinn symmetry 
remains unbroken in the early universe. 



^Work supported by the Department of Energy, contract DE-AC02-76SF00515. 



Contents 



1 Introduction 2 

2 Relic Axions 5 

2.1 Formation of the Axion Condensate 6 

2.2 Small Initial Misalignment Angle 9 

2.3 Dilution of the Axion Condensate 10 

2.3.1 Late Decaying Particle 11 

2.3.2 Saxion Decay 13 

2.4 Inflationary Axion Fluctuations 15 

3 Axion Signatures in the CMBR 17 

3.1 Gravity Waves 18 

3.2 Isocurvature Perturbations 19 

3.3 Non-Gaussianity 25 

4 Constraints on the Axion from the CMBR 30 

5 Remaining Axion Windows 38 

6 Conclusion 41 
A The Model Independent String Axion 43 
B Zero and Finite Temperature Axion Mass 44 
C Parameterizing Non-Gaussian Fluctuations 46 



1 



1 Introduction 



M-theory and perhaps its perturbative string theory hmits represent the best hope for 
a fundamental theory which integrates matter fields, gauge interactions, and gravity. 
This theory makes many interesting predictions for physics at the string or eleven- 
dimensional Planck scales. These scales could in principle be anywhere between 
the electroweak and four- dimensional Planck scales. However, the only compelling 
indirect hint for these scales comes from the unification of gauge couplings due to four- 
dimensional renormalization group running below the unification scale, Mu ~ 2 x 10^^ 
GeV [1]. In this case the fundamental scale is so far above any energy scale available 
in the laboratory that direct probes of string/M-theory seem practically impossible. 
However, energies in the early universe might have been comparable to these scales. 
For example, the energy density during infiation could have been as large as the 
unification scale. It is therefore worth investigating in detail whether any features of 
string/M-theory could give effects in the very early universe which might be detected 
in the present and coming generation of precision cosmological measurements. 

Previous discussions have been rather pessimistic as to the possibility of discern- 
ing effects from such large energy scales with cosmological measurements [2]. How- 
ever, in this paper we show that there is at least one feature of many string and 
M-theory vacua which can be probed by precision cosmological measurements - a 
string/M-theory QCD axion. This conclusion turns out to be nearly independent of 
any cosmological assumptions such as the initial axion misalignment angle, how much 
of the dark matter is relic axions, or dilution of relic axions by a late entropy release. 

A generic feature of supersymmetric string and M-theory compactifications is 
the existence of complex moduli fields. The compact components of these moduli 
ultimately arise in one description or another from p-form gauge fields on non-trivial 
p-cycles in an internal manifold. The vacuum angles for low energy gauge groups 
often transform inhomogeneously under shifts of these compact moduli. In the case of 
QCD such a compact modulus then gives a realization of the Peccei-Quinn symmetry 
in Nambu-Goldstone mode [3]. Supersymmetry breaking generically lifts moduli. 
While not guaranteed, it is possible that the compact components of moduli are 
protected from obtaining a potential from supersymmetry breaking to a fairly high 
order by approximate or exact discrete symmetries [4]. In the case of QCD the leading 
contributions to the potential for such a compact modulus can then arise from infrared 
QCD effects. Such a compact modulus then acts as an invisible axion and provides 
an elegant solution to the strong CP problem. 

In the classes of vacua mentioned above in which gauge coupling unification arises 
from four-dimensional renormalization group running, the Peccei-Quinn scale for a 
string/M-theory QCD axion is generally of order the unification scale, ~ Mu- 
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For example, as presented in appendix A, the Peccei-Quinn scale for the canonically 
normalized model-independent string axion is 

/„ = a/2 ^Mp ~ 10^^ GeV (1) 

4:71 

where au = gfj/'i'n: is the unified value of the fine structure constant and Mp ~ 2.4 x 
10^® GeV is the four-dimensional reduced Planck mass. We therefore focus on Peccei- 
Quinn scales fa ~ 10^^ GeV, but consider the general case of large /„. The interactions 
of such axions are exceedingly weak, first because of the derivative coupling, and 
second because of the large Pcccci-Quinn scale. Even if a relic condensate of such an 
axion comprised all the dark matter it seems practically impossible to detect directly 
in laboratory experiments. 

Precision cosmological measurements do, however, provide the possibility to probe 
a string/M-theory QCD axion. Quantum fiuctuations are imprinted on such an axion 
during inflation. These fluctuations eventually contribute to temperature fluctuations 
in the Cosmic Microwave Background Radiation (CMBR) with the distinctive char- 
acteristics of isocurvature [5] and non-Gaussian [6] components. The magnitude of 
these fiuctuations are controlled by the energy scale of infiation. At present there 
is no direct experimental indication for this scale. However, upcoming experiments 
sensitive to the polarization of the CMBR, such as the PLANCK Surveyor [7] or 
CMBPol could conceivably detect background primordial gravitational waves. The 
presence of a gravitational wave signal in CMBR fiuctuations would represent a true 
triumph for CMBR experiments, and if interpreted as arising from infiation would 
imply a rather high infiationary energy scale. Most importantly for the present dis- 
cussion, it would establish the magnitude of primordial axion fiuctuations imprinted 
during infiation. And for a given Peccei-Quinn scale, /„, this would provide a lower 
limit on the isocurvature and non-Gaussian components in the CMBR fiuctuations. 
So if gravity waves are detected in future CMBR fiuctuations experiments, then the 
existence of a string/M-thcory QCD axion can be tested by a search for isocurvature 
and/or non-Gaussian components also in CMBR fiuctuations. 

The notion that information might be gained about axions from the CMBR is 
not new [6,8-13]. However, these works either assume that the axion makes up all of 
the dark matter (or even the entire critical density) [8,9,11-13], do not allow for the 
possibility of dilution of the axion condensate by a late entropy release [6,8-13], and 
do not consider separately isocurvature and non-Gaussian contributions to CMBR 
temperature fiuctuations [6,8-13]. 

The goal here is to assess how well precision CMBR measurements can probe a 
string/M-theory QCD axion independent of any cosmological assumptions. In order 
to do this the correlation between the magnitude of infiationary gravity wave con- 
tributions to CMBR fiuctuations, and axion induced isocurvature and non-Gaussian 
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components of CMBR fluctuations must be derived for the widest possible range 
of cosmological assumptions: First, the initial axion misalignment angle must be 
taken to be arbitrary within allowed constraints. Second, no assumption must be 
made about how much of the dark matter is made of relic axions. And third, any al- 
lowed, but otherwise arbitrary, late entropy release which could dilute any relic axions 
must be allowed. Obtaining a definitive statement about how well precision CMBR 
measurements can probe a string/M-theory QCD axion then amounts to employing 
the cosmological assumptions which minimize the axion induced isocurvature and 
non-Gaussian contributions to the CMBR fluctuations for a given magnitude of in- 
flationary gravity wave contributions to CMBR fluctuations and Peccei-Quinn scale. 
This corresponds to assuming a vanishingly small average misalignment angle, with 
the relic axion condensate arising only from inflationary induced fluctuations of the 
axion field [10], along with the maximum allowed dilution of this condensate by a late 
entropy release. In this case, any fluctuation induced relic axions turn out to make up 
at most a very small fraction of the dark matter. If this most pessimistic of cosmolog- 
ical assumptions can be probed experimentally, then all possible assumptions can be 
probed: the observational effects of an axion are necessarily larger for non-vanishing 
average initial misalignment angle and non-maximal dilution by an entropy release. 
The motivation for employing the widest possible set of cosmological assumptions 
is that in the era of precision cosmological measurements, ideally it should not be 
necessary to make ad hoc untcstable assumptions about early universe cosmology in 
order to obtain deflnitive tests. 

Even with the above assumptions which minimize the observational effects of relic 
axions we come to a strong conclusion: 

A string/M-theory compact modulus can not play the role of a QCD axion 
and solve the strong CP problem if gravitational waves interpreted as aris- 
ing from inflation with a Hubble constant ifinf ^ lO-*"^ GeV are observed 
by the PLANCK polarimetry experiment. 

This conclusion requires only the current bounds on isocurvature and non-Gaussian 
contributions to CMBR fluctuations and relies on the expectation that the maxi- 
mum possible Peccei-Quinn scale for a string/M-theory axion is the four-dimensional 
Planck scale fa ^ Mp [14]. In terms of an arbitrary Peccei-Quinn scale these bounds 
would require at least /a ^ 3 x 10^° GeV if PLANCK observed gravity waves corre- 
sponding to an inflation scale of i?inf ^ 10^^ GeV, and at least /a ^ 5 x 10^^ GeV if 
gravity waves were observed just below the current WMAP bound of iJi„f < 3 x 10"^^ 
GeV [15]. These bounds do require the mild assumption that the Peccei-Quinn sym- 
metry remains unbroken in the early universe. They would be further strengthened 
by (likely) improvements in the bound on the isocurvature components of CMBR 
fluctuations. 
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In the next section the production of rehc axions is reviewed in all the relevant 
cosmological scenarios including the effects of inflationary induced axion fluctuations 
and the possibility of dilution by an entropy release due a late decaying particle. The 
specific case of the late decay of the saxion superpartner of the axion is considered, 
but shown to generally lead to additional cosmological problems rather than a suc- 
cessful dilution of the axion condensate. The standard results are also extended to 
the case of large Peccei-Quinn scales. In section 3 the isocurvature and non-Gaussian 
contributions to CMBR temperature fluctuations arising from inflationary axion fluc- 
tuations are presented. The current bounds and possibilities for future improve- 
ments on the CMBR probes of the Hubble constant during inflation and isocurvature 
and non- Gaussian components are also discussed. In section 4 the rather stringent 
bounds on the Peccei-Quinn scale which would result from the current isocurvature 
and non-Gaussian CMBR temperature fluctuation bounds are presented in the case 
that primordial gravity waves interpreted as arising from inflation are in fact observed 
in future experiments. The implications for a string/ M-theory QCD axion are also 
considered. In section 5 the scenarios for realizing the Peccei-Quinn mechanism to 
which the stringent bounds would apply are spelled out. Possible loop holes for these 
bounds are also addressed. In section 6 the strong conclusion that a string/M-theory 
QCD axion is not consistent with the observation of primordial gravity waves by the 
PLANCK polarimetry experiment is reviewed. And a corollary is presented regarding 
the implications of observation of isocurvature and/or non-Gaussian components of 
CMBR temperature fluctuations would have if interpreted as arising from a QCD 
axion. 

The relation between the unified fine structure constant and properly normalized 
Peccei-Quinn scale for the model independent string axion is derived in appendix A. 
The parameterization of the finite temperature axion mass resulting from the dilute 
instanton calculation of the ^qco-dependence of the finite temperature free energy 
is given in appendix B. The relation between the dimensionless skewness defined 
in section 3.3 and a parameterization of non-Gaussian components of temperature 
fiuctuations commonly used in the literature is derived in appendix C. 

2 Relic Axions 

Relic axions can be produced in the early universe through a variety of mecha- 
nisms: thermal production, axion string radiation, or condensate formation from 
misalignment of the initial QCD vacuum energy [16, 17, 18, 19]. For the large Peccei- 
Quinn scales relevant here, thermal production is unimportant. In addition, since any 
string/M-theory modulus which contains a QCD axion is present before and during 
infiation any axion strings are infiated away. This leaves coherent production from 
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misalignment as the dominant mechanism for producing rehc string/M-theory QCD 
axions. 

The formation of a QCD axion condensate is reviewed in section 2.1, with special 
attention paid to the case of large Peccei-Quinn scale relevant to a string/M-theory 
axion. It is shown that there are sizeable uncertainties in the calculation of the axion 
relic density due to strong QCD effects for Peccei-Quinn scales in the region of interest, 
but that this does not affect the final conclusions. An updated version of the classic 
cosmological bound on the Peccei-Quinn scale for an average misalignment angle is 
presented taking into account recent determinations of the dark matter density. This 
bound is significantly exceeded for fa/N ~ 10^^ GeV, the value that might be expected 
for a string/M-theory QCD axion. So the possibility of accommodating such an axion 
by reducing relic production with a small misalignment angle is considered in section 
2.2. The complimentary or additional possibility of reducing relic axion density by 
a late entropy release is discussed in section 2.3. It is shown that the maximum 
possible dilution by a late decaying particle is alone not sufficient to accommodate 
fa/N ~ 10^^ GeV with average misalignment angle. Such an axion requires a small 
misalignment angle. The inflationary quantum fluctuations imprinted on an axion 
are reviewed in section 2.4. These fluctuations are shown to lead to a minimum relic 
density of axions which depends on the Hubble constant during inflation and the 
Peccei-Quinn scale. 



2.1 Formation of the Axion Condensate 

A relic axion condensate is formed in the early universe from misalignment of the 
initial QCD vacuum angle [16] when the axion begins to oscillate at ?>H ~ ma(Tosc), 
where H is the Hubble parameter. Given the temperature dependence of the axion 
mass ma{T), this condition along with the Friedman equation, 3H'^Mp = p, and 
thermal energy density in a radiation dominated era p — [n'^ / 30) g^T^^^ deflnes the 
temperature, Tosc, at which the axion begins to oscillate 

Tosc ^ (10/ gyf\ma{Tosc)Mpf^ (2) 

For oscillation temperatures just above Aqcd, 9* — 61.75. The parameterization of 
the high temperature dependence of the axion mass for T ^ Aqcd is reviewed in 
Appendix B, and may be written [20, 17] 

^(T) ^ ^ c ( J^SS^\" ( for T > Aqcd (3) 

' rua V200 MeV; \ T ) ~ Qod \ ) 
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where C ^ 0.018 is a parameter defined in Appendix B. Witfi tfiis, the oscillation 
temperature is then 

This expression for the oscillation temperature is only appropriate for Tqsc ^ Aqcd 
corresponding to fa/N < 2 x 10^^ GeV. The case of T^^c ^ Aqcd for larger Peccei- 
Quinn scales is considered separately below. 

The local axion number density in the condensate when it is formed is 

ua ^ u\i{T..c)ma{fa/Nm + al)f{9^) (5) 

where = '^.^(O) is the zero temperature mass, /c is a correction for temperature 
dependence of the axion mass, /c — 1.44 [17] for the ^(T) parameterization (3), and 
f{x) is a correction for anharmonic effects of the axion potential; for x -C 1, /(x) — > 1. 
And 9i = {9) is the initial QCD vacuum angle zero-mode averaged over the presently 
observable universe, and (Tq = {{6 — {9))"^) is the mean square fluctuations of the 
initial QCD vacuum angle, produced for example during an early epoch of inflation. 
Since the number density redshifts like the entropy density, the local axion density 
today is given by 

Pa - {rial s)maSQ-i (6) 

where s — {27:'^ /A5)g^sTosc thermal entropy density when the axion begins to 

oscillate, sq — (27r^/45)g'*so7'o thermal entropy density today with g^so — 3.91 

and To ~ 2.73 K. The factor 7 is a possible dilution factor due to entropy release 
after the axion begins to oscillate, 7 = Sosc/Sq- Using all this with the oscillation 
temperature (4) leads to a relic axion density, fla = Pa/ {'^H^Mp), of 

where h = Hq/(100 km s^^ Mpc~^). This expression for the relic axion density is 
only valid for at least fa/N < 2 x lO^^GeV. 

Since the expression (7) for the relic axion density uses the high temperature 
expression (4) for the oscillation temperature, it does not apply for T^sc ^ Aqcd- For 
oscillation temperatures in this range (corresponding to large Peccei-Quinn scales), 
the axion begins to oscillate in effectively the zero temperature potential and another 
expression for the relic density is obtained. For oscillation temperatures just below 
Aqcd, 9* — 10.75. With this, and for a general temperature dependence for the axion 
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mass, ^{T), the oscillation temperature (2) is 

/iq16 r<pv\ 

7;,e~ 950 MeV ( j j ^(Tosc)'/' for To^e < Aqcd (8) 

This expression for the oscillation temperature is only appropriate for T^sc ^ ^qcb, 
corresponding to at least fa/N > 2 x 10^''' ^(Tqsc) GeV. For Tqsc ^ Aqcd, the rehc 
axion density from (5) and (6) is 

^ah' c 5 X 10^ [ jil^Gev ] (^^' + ^^') ^^^^^^ MT^^'^y (9) 

While there is some residual temperature dependence of the axion mass for T^sc ^ 
Aqcd, this is unimportant and fcC{Tosc)~^ ~ 1 is a good approximation well into this 
regime. The relic density (9) is only valid for at least fa/N > 2 x 10"'^'^ GeV. Note 
that is has a different parametric dependence on fa/N than (7) appropriate for small 
Peccei-Quinn scales. 

Note that the relic axion densities (7) and (9) calculated for oscillations tempera- 
tures (4) and (8) above and below Aqcd do not have overlapping regions of validity. 
This indicates that in the transition region, 2 x 10^^ GeV < fa/N < 2 x 10^'^ GeV, 
neither expression is completely accurate. This is in fact to be expected for Tqsc ~ 
Aqcd- The high temperature expression (4) is derived using a dilute instanton gas 
approximation [20] as described in Appendix B. This likely receives non-trivial cor- 
rections from multi-instantons for Tqsc somewhat above Aqcd- Likewise the low tem- 
perature expression (8) likely receives non-trivial corrections from spontaneous chiral 
symmetry breaking effects for Tqsc just slightly below Aqcd- It is perhaps unfortunate 
that for a Peccei-Quinn scale relevant to a string/M-theory axion, strong QCD effects 
are important during formation of the axion condensate. However, for fa/N well out- 
side the transition region the true expression for the axion relic density asymptotes 
to the expressions given above, and should smoothly interpolate between these lim- 
iting forms. And the bounds described below turn out to apply to all Pcccci-Quinn 
scales up to a very large value which is well outside the transition region. So the 
final conclusions are not sensitive to unknown strong QCD effects. The relic densities 
(7) and (9) turn out to be equal for fa/N ~ 6 x 10^^ GeV. This is in fact outside 
the transition region given above, due mainly to the abrupt change in the number of 
degrees of freedom at a temperature T ~ Aqcd- In what follows we simply use (7) 
for Peccei-Quinn scales lower than this value and (9) for higher scales. Also note that 
the expressions for the relic densities (7) and (9) only apply if the universe is radia- 
tion dominated when the axion condensate forms. The case of matter domination is 
presented in section 2.3. 

Most of the relic axions have momenta much less than the zero temperature axion 
mass. Pa <S rua, and so act as cold dark matter (CDM). The present CDM energy den- 
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sity has been accurately determined in a ACDM cosmology by recent measurements 
to be f2cDM^^ = 0.113 ± 0.010 [21]. At most, relic axions could saturate this density 
and make up all the CDM. Alternately, the CDM might consist of more than one 
species, with relic axions contributing only part of the total observed CDM density. 
So the relic axion density is bounded from above by 

^ah^ < 0.12 (10) 

The average mean square initial misalignment angle over all inflationary domains 
is {9'^) = 7r^/3, for which the anharmonic correction is /{tt'^/S) ~ 1.2 [17]. So for the 
average misalignment angle the relic axion density (7) without any dilution, 7 = 1, 
along with the bound (10), require that the Peccei-Quinn scale is bounded from above 

by 

fa/N < 3 X 10^^ GeV for {9^) = n^S and 7 = 1 (11) 

This is somewhat smaller than the canonical fa/N < 10^^ GeV bound arising from 
flah"^ ^ 1 quoted in the older literature [16] mainly due to the recent improvement 
in measurement of the CDM density quoted above. 

2.2 Small Initial Misalignment Angle 

A QCD string/ M-theory axion with fa/N ~ 10^^ GeV is clearly inconsistent with the 
bound (11) for the average misalignment angle and without any dilution of the axion 

condensate after its formation. One possibility which could accommodate such an 
axion is that the average misalignment angle within the observable universe is small. 
For the relic axion density (7) with 7=1, the bound (10) then requires 

{9^ + a^,y/^ < 2x10-^ (12) 

In fact, with random initial conditions prior to an early epoch of inflation, there are at 
present distant regions of the universe which sample all values of initial misalignment 
angles. So it is possible in principle that 9f + aQ happens to be sufficiently small in 

our observable universe to allow fa/N ~ 10^^ GeV. Regions with misalignment angles 
which are larger (smaller) than ours then have a larger (smaller) Hubble constant at a 
given temperature. 

In this scenario the question naturally arises though as to why our initial mis- 
alignment angle is so small. Linde has argued that with as large as considered 
here there may be a dynamical selection effect which makes the observation of small 
misalignment angle likely [9] (for a recent discussion see [22]). At fixed large /„, the 
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density of gravitationally collapsed objects such as galaxies grows as a rapid power of 
6f [9]. Since very dense galaxies may be inhospitable to observers such as ourselves, it 
may then be more likely that a given observer measures a small misalignment angle. 
However, quantifying the probability distribution of observers as a function of galac- 
tic density remains an open problem. In addition, such arguments could not explain 
an initial misalignment angle significantly smaller than the bound (12). 

Independent of the question of possible selection effects for the observed value 
of the initial misalignment angle, it is possible in principle that the misalignment 
angle in the observable universe is so small for whatever reason that relic axions 
comprise a vanishingly small fraction of the CDM. In this case it might appear that 
there are no observational effects. However, quantum fluctuations of the axion during 
inflation lead to a non- vanishing mean square fluctuation of the QCD vacuum angle, 
ae 7^ 0, as discussed in section 2.4. This gives a lower limit on the relic axion 
density (7) or (9), and therefore on the observable effects. The magnitude of these 
effects depends crucially on the Hubble constant during inflation, as discussed below. 
The goal here is to estabhsh how well a QCD string/ M-theory axion can be probed 
by precision cosmological measurements independent of any assumptions about the 
initial misalignment angle. To do this, the Hubble constant during inflation must 
therefore be measured by independent means in order to establish a minimum relic 
axion density, as discussed in section 3.1. 

2.3 Dilution of the Axion Condensate 

Another possibility which might help accommodate a QCD string/M-theory axion 
is that there is a late entropy release with a reheat temperature Trh < ^osc which 
dilutes the axion condensate after its formation [23]. For the relic axion density (7) 
with average misahgnment angle {9'^) — 7r^/3, the rehc density bound (10) requires a 
dilution factor of 



A late entropy release must satisfy a number of non-trivial constraints in order 
to successfully dilute the axion condensate. First, the entropy release must occur 
after the condensate is formed at Tosc ~ Aqcd- An entropy release at a higher 
temperature does not modify the relic axion density. Second, the entropy release 
must not adversely modify the successful predictions of relic light element abundances 
from big bang nucleosynthesis (BBN). The ^Hc abundance is determined in standard 
radiation dominated BBN mainly by the primordial baryon to entropy ratio, rib/s, and 
the neutron to proton ration, n/p, when the weak interactions drop out of equilibrium 
at T ~ 1 MeV. Standard BBN then generally requires that any entropy release 




(13) 
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terminate in a radiation domination era with a reheat temperature of at least Trh ^ 
1 MeV and possibly higher. Finally, the entropy release must not over-dilute the 
primordial relic baryon density. 

The case of an entropy release by a late decaying particle is considered below. 
The special case of a late decaying saxion superpartner of the axion is also presented. 
The very unlikely possibilities of a first order phase transition or very late inflation 
which could dilute the axion condensate are commented on in section 5. 

2.3.1 Late Decaying Pcirticle 

A very natural way to obtain a late entropy release is by the late decay of a massive 
non-relativistic particle. In order to have a non-trivial dilution of the axion condensate 
the late decaying particle must dominate the energy density for at least Tqsc ^ Aqcd- 
This can occur in a number of ways. For a scalar particle, a relic condensate formed 
from an initial misalignment of the scalar field from the minimum of the potential 
can dominate at a very early epoch. For either a scalar or fermion particle, relic 
particles produced during reheating after inflation, or simply from thermal freeze out 
can dominate at an early epoch. For the discussion here, detailed models for the 
origin of the late decaying relic particles are not necessary. The dilution of the axion 
condensate can be parameterized solely in terms of the particle reheat temperature 
as described below. 

For a late decaying particle the BBN bounds on the reheat temperature are shghtly 
more stringent than the case of a general entropy release discussed above. An ex- 
ponentially small number of late decaying particles survive to temperatures below 
the reheat temperature. If hadronic decay channels are available (which is inevitable 
in the thermalization process) then these very late decays can enrich n/p even for 
^RH ^ 1 MeV. And this can be a signiflcant effect since n^/s is so small. For a 
late decaying particle with m ~ 1 TeV, a reheat temperature of Trh ^ 6 MeV is 
required in order to ensure that the "^He abundance is not too large [24]. So in order 
to dilute the axion condensate while keeping the successful predictions of BBN a late 
decaying particle must flrst, dominate the energy density for at least T > Aqcd and 
second, have a reheat temperature in the range 6 MeV ^ Trh ^ Aqcd- 

The rehc axion density which results with the late decay of a massive particle 
can be calculated from the axion fractional energy density when the axion condensate 
is formed. If a massive non-relativistic particle dominates the energy density, the 
universe is matter dominated at this epoch with = n(pm(p ~ SH^M^. The number 
density in the axion condensate when it forms at 3H ~ ma(Tosc) is given by (5). 
Since the axion and massive particle number densities both redshift like matter, the 
ratio Ua/n^ is constant after the axion condensate is formed and before the massive 
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particle decays. Using all this, the ratio of axion density pa — maUa to massive 
particle density during this epoch is 

^ - /e|e(Tosc)-^^^ + m) (14) 

When the massive particle decays its energy is converted to thermal radiation. In 
the sudden decay approximation the resulting thermal energy and entropy densities 
are ~ {n'^ /30)g^T^^ and s ^ (27r^/45)5f*T^jj. The relic axion number density to 
entropy density ratio at this epoch is then 

ria _ pa/mg 
s P4> s 

with p^j s (3/4)Trh- Since the number density redshifts like the entropy density 
after the late decay, the local axion density today is given by (6) using the ratio (15). 
Combining all this, the relic axion density with the late decay of a massive particle 
is [19] 

Note that in this scenario appropriate to a matter dominated era, the relic axion 
density has a parametric dependence on fal^ which is different than both (7) and 
(9) appropriate to a radiation dominated era. 

Starting from a very early radiation dominated era, the temperature for a given 
Hubble parameter during the subsequent late decaying particle matter dominated era 
is smaller than that for the same Hubble parameter of a radiation dominated era. So 
as long as the late decaying particle dominates the energy density at a temperature 
at least somewhat above Aqcd, the temperature that the axion begins to oscillate, 
Tosc, is at least somewhat below Aqcd- The axion condensate then forms in the zero 
temperature potential, and /c'C(^osc)^^ ~ 1 is a good approximation. This implies that 
if the late decaying particle comes to dominate the energy density for a temperature 
well above Aqcd (which is likely in most scenarios) then the relic abundance (16) does 
not suffer unknown uncertainties from strong QCD effects in the region of interest, 
unlike the relic densities (7) and (9) appropriate to a radiation dominated era. 

For the average mean square misahgnment angle of (^^) = 7r^/3, the rehc density 
(16) with maximum dilution, Trh — 6 MeV, along with bound (10), requires that the 
Peccei-Quinn scale is bounded from above by 

fjN ^ 3 X 10^^ GeV for If) = ttVS and Maximum Dilution (17) 

A string/M-theory QCD axion with /a/A^ ~ 10^*^ GeV is inconsistent even with this 
maximum dilution bound. So dilution from late decay of a massive particle turns 
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out alone to be insufficient to accommodate such an axion with generic misahgnment 
angle. A small misalignment angle is therefore required in any cosmological scenario 
with a Peccei-Quinn scale which far exceeds the bound (17). 



2.3.2 Saxion Decay 

In any supersymmetric model the pseudoscalar axion has a scalar saxion superpartner. 
For the model independent string axion the saxion is just the dilaton. In the super- 
symmetric limit, Peccci-Quinn symmetry ensures that the entire axion supermultiplet 
remains masslcss to all orders in perturbation theory. However, with supersymmetry 
and U{l)ji breaking, only the axion mass is protected by Peccei-Quinn symmetry, 
and the saxion potential is determined by supersymmetry breaking. The saxion mass 
is then naturally much larger than the axion mass. A coherent condensate of the 
saxion can form at a very early epoch from saxion misalignment and dominate the 
energy well before the axion condensate forms. The saxion condensate eventually de- 
cays, with the decay products rethermalizing to a standard radiation dominated era. 
So a supersymmetric axion model naturally contains all the elements which could in 
principle dilute the axion condensate [25, 26]. However, as shown below, the reheat 
temperature is generally too low for fa/N ~ 10^^ GeV. This complicates rather than 
improves the cosmological scenario. 

The detailed cosmological scenario which is obtained with a supersymmetric axion 
model depends on the saxion reheat temperature. This is determined by the saxion 
decay width, F, and mass, and is therefore sensitive to details of supersymmetry 
breaking. In the sudden decay approximation the saxion condensate decays when 
F ~ if with its energy converted to thermal radiation, (7r^/30)fir^T^^ ~ 3F^Mp. The 
decay width depends on the saxion couplings which are model dependent. However, 
if the anomalous axion-gluon coupling (74) given in Appendix A arises directly in 
the high energy theory (rather than from integrating out low energy standard model 
fields) as is the case for large classes of string/M-theory axions including the model 
independent string axion, the saxion-gluon coupling is determined by supersymmetry 

where the gluon fields are canonically normalized. This coupling gives a decay rate 
for s ^ gg oi 

oP' rrfi 

In the sudden decay approximation the decay rate (19) gives a reheat temperature of 

. 1 X 10- MeV (^) (20) 
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with n = 2 corresponding to the two derivatives in the couphng (19). As another 
example, the DFS [27] class of axion models couple only to the visible sector through 
the Higgs sector. The couplings to standard model fields then come from mixings with 
Higgs bosons. In this case the leading coupling of the saxion are to heavy quarks 

ruQQQ (21) 



fa/N 



For the DFS axion c = | for each quark and = 6. Couplings to gluons include an 
additional loop factor and are therefore suppressed. The coupling (21) gives a decay 
rate for s — > of 




r(^-««-i6^a7i#i'--^i '''' 

For a DFS saxion heavier than the top quark, the dominant decay is to tops, s ti, 
with a reheat temperature from (22) which numerically turns out to be essentially 
identical to (20) with n — corresponding to the absence of derivatives in the coupling 
(21). 

For a saxion mass of ~ 1 TeV, the reheat temperature (20) is consistent with 
the standard BBN requirement of Trh ^ 6 MeV [24] only for at least fa/N < 
2 X 10^^ GeV. So for fa/N ~ 10^^ GeV appropriate to a string/M-theory QCD axion, 
the rather late decay of a saxion condensate formed at an early epoch leads to an 
additional saxion problem which must be solved, rather than contributing successfully 
to dilution of the axion condensate. A very large saxion mass could however increase 
the reheat temperature sufficiently to be consistent with BBN. For example, in the 
case of an saxion with high energy coupling (18) and Peccei-Quinn scale fa/N ~ 
10^^ GeV, a saxion mass of iris ^ 70 TeV would just barely give Trh ^ 6 MeV. 
Although as discussed in section 2.3.1, dilution by late decay alone is not sufficient to 
accommodate such an axion without small initial misalignment angle. And obtaining 
such a heavy saxion seems difficult in most models of supersymmetry breaking with 
superpartners just above the electroweak scale. In many models of supersymmetry 
breaking the saxion mass is in fact much lighter than assumed above, <^ TeV, 
leading to an even more severe saxion problem. In addition, the axino fermionic 
partner of the axion has, for a given mass, a similar decay rate and reheat temperature 
to those given above [28]. So if sufficiently produced, relic axinos can also lead to a 
problem. However, in some scenarios for supersymmetry breaking the axino can be 
the lighest super symmetric particle, and might be a CDM candidate [28]. 

Any complete cosmological scenario which includes a supersymmetric QCD axion 
must address the potential problems associated with the saxion and axino partners, 
presumably by limiting the initial cosmological production of these states. Since 
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the goal here is to assess how well precision cosmological measurements can probe 
a string/M-theory QCD axion independent of specific cosmological assumptions, de- 
tailed models which address the saxion and axino problems will not be presented. 
Instead, we show below that in any model in which these problems are in fact solved, 
precision CMBR measurements can in principle definitively probe such an axion. 

2.4 Inflationary Axion Fluctuations 

The relic axion density is proportional to the quadrature sum of the zero-mode and 
root mean square fiuctuations of the initial QCD vacuum angle, 6f + aj. This combi- 
nation amounts to an arbitrary initial condition for the mean square axion field when 
the axion condensate is formed. But in an infiationary cosmology the mean square 
axion fluctuations are related to the Hubble constant during inflation. In this sec- 
tion, inflationary production of axion fluctuations are reviewed. The minimum relic 
axion densities associated with these fluctuations [10] for the scenarios described in 
the previous sections are presented. These minimum relic axion densities are crucial 
in determining a lower limit on the observational effects described in the subsequent 
sections. 

All massless fields undergo de Sitter fiuctuations during an infiationary epoch. 
At zero temperature the leading Peccei-Quinn breaking effects which lift the axion 
potential must be due to QCD instantons in order to obtain a successful solution of the 
strong CP problem. These effects are exponentially suppressed during infiation; and 
as long as the Peccei-Quinn symmetry is not explicitly broken in the early universe, 
the axion is effectively massless during inflation. The de Sitter power spectral density 
of fluctuations of a massless fleld such as the axion 



where the k-space density is VcPk/{2TT)^ — Vdkk"^ /27r'^ and V is the quantization 
volume. The power spectral density (23) is nothing other than the mean square 
fiuctuations of the field, Va = cr"^- So the root mean square fiuctuations of the axion 
field and initial misalignment angle are related to the Hubble constant by 




(23) 



is related to the Hubble constant during inflation by 




(24) 



(Je = 



Mfa/N) 



(25) 



15 



where Tqh is the Gibbons-Hawking temperature [29]. Numerically, the variance of 
the initial misalignment angle is related to the Hubble constant during inflation and 
Peccei-Quinn scale by 

ae ~ 1.6 X 10 M ,,,, -^7^ (26) 



1013 GeVj [ fjN I 

This inflationary produced variance provides a lower bound on the combination Of+UQ 
that appears in the expressions for the axion relic density. 

In order to establish a lower limit for the observational effects of relic axion fluc- 
tuations in the CMBR discussed in the next section, it is useful to determine the min- 
imum relic axion density as a function of Hubble scale during inflation for each late 
time cosmological scenario described in the previous subsections. Consider first the 
case in which the axion condensate forms during a radiation dominated era without 
any dilution by a late decaying particle. For a Peccei-Quinn scale of /„/ A?" < 6 x 10^^ 
GeV, the relic density (7) with C = 0.018 and Aqcd = 200 MeV, and for an aver- 
age initial misalignment angle which is small compared with the variance induced by 
inflationary fluctuations (26), namely 9f < aj, is bounded by 

With the same late time cosmological scenario, but with at least fa/N ^ 6 x 10^^ 
GeV, the relic density (9) with /c'C(^osc)~^ — 1 and with small average initial mis- 
alignment angle, 9f ^ a^, gives a bound on the relic axion density of 

The lower bounds (27) and (28) for small and large Peccei-Quinn scales respectively 
depend on inverse powers of fa/N because the variance in initial misalignment angle 
(25) is inversely proportional to the Peccei-Quinn scale. And both scale hke the 
square of the Hubble constant during inflation since the relic density scales like the 
square variance of axion fleld, and the variance scales like the Hubble constant. 

Next consider the case of a late entropy release in which the axion condensate is 
formed during a matter dominated era with a later transition to standard radiation 
domination by reheating from a massive particle. In this case the relic density (16) 
for an average initial misalignment angle which is small compared with the variance 
induced by inflationary fluctuations (26), namely 9^ < ag, is bounded by 
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Note that this bound on the residual rehc axion density induced by inflationary 
fluctuations is independent of the Peccei-Quinn scale because in this case the relic 
density (16) is proportional to [fa/NY while the variance in initial misalignment 
angle (25) squared is proportional to {fa/N)~'^. 

The bound (29) is only applicable if the axion condensate is formed during the mat- 
ter dominated era before the dominating particle decays, which requires 3if (Trh) ^ 
ma(7RH)- For a small enough axion mass the axion condensate is formed in the 
radiation dominated era subsequent to the late particle decay. This will occur for 
3H{T^ji) > ma{T^B.) corresponding to a Peccei-Quinn scale 



In this extreme case the relic density (9) for formation of the axion condensate in 
the zero temperature potential in a radiation dominated era is applicable. And the 
bound (28) on the relic axion density from inflationary induced axion fluctuations 
is obtained. Such a large Peccei-Quinn scale would seem unnatural and does not 
occur in any known string/M-theory compactifications [14]. However, it is useful to 
consider this regime in order to assess the strength of the bounds described in section 
4. 



A string/M-theory QCD axion with fa/N ~ 10"*^^ GeV is likely to be practically 
impossible to detect in the laboratory. However, inflationary produced fluctuations 
of such an axion do give rise to two distinctive signatures in CMBR temperature 
fluctuations which are in principle measurable. The flrst is an isocurvature component 
of temperature fluctuations [5] which are discussed in section 3.2. The second is a 
non-Gaussian component of temperature fluctuations [6] discussed in section 3.3. 

Both of these signatures depend on the relic axion density. The minimum residual 
relic densities arising from inflationary axion fluctuations presented in section 2.4 are 
therefore a crucial ingredient in establishing lower bounds for the observable isocur- 
vature and non-Gaussian effects in the CMBR temperature fluctuations. These lower 
bounds, however, depend on the Hubble constant during inflation. So it is necessary 
to have some independent measurement which can in principle establish the Hubble 
constant during inflation. This can in fact be accomplished by a measurement of pri- 
mordial inflationary gravity wave contributions to CMBR temperature fluctuations 
as described section 3.1 [30]. 




(30) 



3 Axion Signatures in the CMBR 
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3.1 Gravity Waves 



All massless fields undergo quantum de Sitter fluctuations during inflation, as dis- 
cussed in section 2.4. This includes the metric as well as an axion. So a measurement 
of the magnitude of metric fluctuations generated during inflation would establish 
the magnitude of the canonically normalized axion fluctuations generated during in- 
flation. And for a given Peccei-Quinn scale this would then establish a lower limit 
for the magnitude of observable effects of a string/M-theory QCD axion in CMBR 
temperature fluctuations described in the following subsections. 

The power spectral density of inflationary dimensionless metric fluctuations may 
be written 

2(|/.+kP + |/^xkP) (31) 

where h+ and hx are the two polarization modes of the metric perturbations. This 
power is related to the Hubble constant during inflation by [15], 

This is identical to the power spectral density of axion fluctuations (24) up to a 
normalization related to canonical normalization of metric fluctuations and overall 
normalization chosen in (31). 

Metric perturbations generated during inflation appear as gravity waves re-entering 
the horizon after inflation. These gravity waves induce tensor B-modes in CMBR tem- 
perature fluctuations. In order to characterize the magnitude of the power in gravity 
waves it is useful to define the ratio of tensor to scalar metric perturbations. Here we 
adopt the normalization used by the WMAP collaboration [15, 31] in which the tensor 
to scalar ratio is defined by the ratio of spectral power densities at some reference 
scale ko, 

r = ^ (33) 

where 

Vn{k) ^ (34) 

is the power spectral density in scalar metric perturbations and TZ is the gauge invari- 
ant dimensionless curvature perturbation [32]. It is worth noting that the definition 
(33) has the advantage that the relation between r and the Hubble constant during 
inflation does not depend on cosmological parameters such as the present dark matter 
or dark energy fractions (in contrast to some other definitions of the tensor to scalar 
ratio in the literature). The scalar curvature perturbation spectral power density 
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measured from CMBR temperature fluctuations by WMAP under the assumption 
that they arise predominantly from scalar metric fluctuations is [15] 



Vn{ko) ^ 2.1 X 10"^ at ko = 0.002 Mpc 



-1 



(35) 



With this, the variance of inflationary induced fluctuations in the initial QCD mis- 
alignment angle (26) may be related to the scalar to tensor ratio (33) using (32) 



For a given Pcccei-Quinn scale this gives a relation between an observable quantity 
and the variance of the initial misalignment angle. 

The current WMAP bound on the scalar to tensor ratio from the non-observation 
of B-mode tensor components in CMBR temperature fluctuations is r < 1.3 [15]. 
With the deflnition (33) and (32) this corresponds to an upper limit on the Hub- 
ble constant during inflation of iJinf ^ 2.8 x 10^^ GeV. Future experiments will 
be capable of probing signiflcantly smaller values of the scalar to tensor ratio r and 
concomitantly smaller values of the Hubble constant during inflation. The PLANCK 
polarimctry experiment could ultimately be sensitive at the one sigma level to at best 
roughly r > 5 x 10~^ corresponding to ifinf ^ 5 x 10^^ GeV [33, 34]. This represents 
a very optimistic sensitivity for a positive observation. So throughout we will consider 
-f^inf ^ 10^^ GeV as a reasonable benchmark for a positive observation of primor- 
dial gravity waves by the PLANCK polarimetry experiment. If achieved this would 
represent a significant improvement over the current WMAP bound. If PLANCK 
does indeed observe tensor B-modes at any level between the current WMAP bound 
and its ultimate sensitivity, this would imply a lower bound on the axion relic den- 
sity for a given Peccei-Quinn scale in the various cosmological scenarios discussed in 
section 2.4. The ultimate sensitivity to primordial gravity waves achievable in fu- 
ture CMBR experiments is limited not only by issues of instrumentation, but also 
from the fact that tensor B-mode fluctuations can be generated at second order in 
primordial scalar fluctuations from gravitational lensing effects. This provides an 
irreducible background for any future probe of primordial gravity waves in CMBR 
experiments. In terms of the scalar to tensor ratio this has been estimated to be 
roughly r ^ 7 x 10"^ corresponding to i^inf ^ 2 x 10^^ GeV [33, 35]. 

3.2 Isocurvature Perturbations 

Quantum fluctuations generated during inflationary can lead to observable signatures 
in CMBR temperature fluctuations in a number of ways. Fluctuations of the inflaton 
field give rise to fiuctuations in the scalar curvature and lead to adiabatic fiuctuations 



by 




(36) 
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in the CMBR. Fluctuations of the metric lead to tensor B-model fluctuations in the 
CMBR, as discussed above. Inflationary induced fluctuations of an axion fleld turn 
out to lead isocurvature fluctuations in the CMBR [5] as reviewed below. This dis- 
tinctive signature provides a probe for the existence of an axion field during infiation. 

Isocurvature perturbations correspond to fluctuations in the local equation of state 
of some species, 5{ni/s) ^ 0, with no fluctuation in the total energy density, 5p — 
These are in some sense orthogonal to adiabatic perturbations which correspond to 
fluctuations in the total energy density, 5p 7^ 0, with no fluctuation in the local 
equation of state, 5{ni/s) = 0. During inflation the potential for a QCD axion is 
essentially flat as long as the Peccei-Quinn symmetry remains unbroken during this 
epoch, as discussed above. Fluctuations in the axion fleld therefore do not affect the 
total density and are isocurvature. 

In order to discuss both adiabatic and isocurvature fluctuations it is useful to 
deflne the fractional fluctuation in the number density of the i-th species divided by 
the entropy density 

5..^ = ^-3^ (37) 

iii/s rii T 

where the second equality follows from s oc . Throughout, all fields other than 
the axion are assumed to undergo adiabatic fiuctuations, = for i 7^ a. Relaxing 
this assumption would only allow for a smaller axion contribution for a fixed total 
isocurvature contribution to CMBR temperature fiuctuations, and so could only result 
in stronger bounds than those given in section 4. As described above, isocurvature 
fiuctuations of the axion on length scales larger than the horizon do not modify the 
total density on these scales, Sp^so — 0. After the axion condensate forms, the total 
energy density is a sum over non-relativistic species including the axion plus radiation 

p = ^ mini + maUa + Pr (38) 

i 

where p^. oc is the radiation energy density, and throughout X]j is over all non- 
relativistic species but the axion. The number density and temperature fiuctuations 
for an isocurvature fiuctuation on scales large compared with the horizon in this epoch 
are then related by 

Spiso = X] l^i^l^i + IT^aSria + Apr — — ^ (39) 

i 

Since all species other than the axion are assumed to fiuctuate adiabatically. Si = 0, 
the definition (37) relates the number density fiuctuation in each species to the initial 
isocurvature temperature fiuctuation by Srii = 3ni{6T/T). With this, the isocurvature 
constraint (39) then relates the isocurvature temperature fiuctuation to the axion 
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number density fluctuations after the axion condensate is formed 



where 5pa = rriaSna. The relation (37) may then be used to write the axion induced 
temperature fluctuation on super-horizon scales in terms of the fractional axion fluc- 
tuation Sa 

^ ~ S (41) 

T S{J:^Pi + Pa)+^Pr " ^ ^ 

This relation is more useful than (40) for arbitrary times after the axion condensate 
forms since Sa rather than Sua/na is approximately constant for fluctuations on super- 
horizon size scales [36]. 

The fractional axion fluctuation Sa can be related directly to the inflationary 
induced fluctuations in the axion fleld. When the axion condensate is formed at 
a temperature of Tqsc ~ Aqcd, in any of the cosmological scenarios discussed in 
section 2, it comprises a small fraction of the total energy density, pa -C J2iPi + Pr- 
From (40) it follows that just after the axion condensate is formed, the axion induced 
fractional temperature fluctuations are much smaller than the fractional axion number 
density fluctuations, {ST/T)init ^ ('5^a/^a)init- Prom this it follows that the fractional 
fluctuation Sa deflned in (37) just after the axion condensate is formed may then be 
related to fluctuations in the initial misalignment angle 

Sua {{e) + Sef - {9') ^ 29,59 + {5ef - al 

" " ua init " {O-") {9^) 9} + al ^^"^ 

where 59 and gq are understood to be due to the initial inflationary induced fluc- 
tuations in the misalignment angle, 9i is the average misalignment angle over the 
currently observable universe, and Ua oc ignoring anharmonic effects (this is a very 
good assumption for /a/A^ ~ 10^^ GeV since ^ <C tt is required in this case). And 
since Sa is approximately constant for super-horizon size fluctuations [36], the rela- 
tion (42) between Sa and the initial misalignment angle fluctuations remains good for 
super-horizon modes for arbitrary times after the axion condensate forms. For later 
reference, with the relation (42), the mean square of the fractional axion fluctuation 
on super-horizon scales is related to the mean square fluctuations of, and average 
value of, the initial misalignment angle by 

(53 - ^ol^t^^ (43) 

where 59 = 9 — (9) is assumed to be Gaussian distributed, 9, = (9), and ag ~ 
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As discussed below, the most relevant axion induced temperature fluctuations are 
those on the largest angular scales. These temperature fluctuations enter the horizon 
well into the matter dominated era for which J2i Pi + Pa Pr- So from (41) the axion 
induced initial temperature fluctuations on large angular scales at horizon crossing 
are related to the axion isocurvature fluctuation Sa by 

T - -It'' 

1 horizon O i 

where Qa = Pa/ Pc is the fractional axion density and Qm = (Si Pi + Pa) / Pc is the 
fractional density of all non-relativistic matter including baryons, axions, and any 
other species which contribute to cold dark matter, = ^cbm + ^b- Ftom WMAP 
the total matter density in a ACDM cosmology is measured to be [21] 

0.135^^09 (45) 

An experimental observation of temperature fluctuations corresponds to the initial 
fluctuation at horizon crossing (44) plus an additional Sachs- Wolfe contribution [37] 
from the integrated redshift between the surface of last scattering and detector. The 
latter contribution turns out to be five times smaller [38] than that from the initial 
temperature fluctuation (44), 

^ ^--^S. (46) 

T sw 15 ilm 

The total observable isocurvature temperature fluctuation is then related to the axion 
isocurvature fluctuation by 

The expression (47) relating axion induced temperature fluctuations with the frac- 
tional axion fluctuation applies strictly only on the largest angular scales. The spec- 
trum of temperature fluctuations is modifled substantially on smaller scales corre- 
sponding to multipoles well beyond the Sachs- Wolfe plateau. This is illustrated in 
Fig. 1 which gives the CMBR temperature fluctuation power spectrum for pure adi- 
abatic or pure isocurvature initial fluctuations in a ACDM cosmology generated by 
CMBFAST version 4.3 [39]. The power in isocurvature fluctuations falls only by a 
factor of two for roughly ^ ~ 50 but continues to fall rapidly for larger £. This rapid de- 
crease in isocurvature contributions to the power in CMBR temperature fluctuations 
compared with adiabatic contributions arises for two reasons. First, the growth of 
sub- horizon size isocurvature and adiabatic temperature fluctuations are completely 
out of phase with one another. So while the magnitude of adiabatic fluctuations 
are growing towards the first peak, the magnitude of isocurvature fiuctuations are 
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Figure 1: Temperature fluctuation power spectrum with pure adiabatic (solid line) or 
pure isocurvature (dashed line) initial fluctuations, AT{9, (f))/T = J2em(^emyim{0, (j)), 
d = {\aim\'^). ACDM cosmology with Qa = 0.73, Qcdm = 0.22, = 0.05, H = 
72 km Mpc~i, and spectral index n—1. 

falling towards a first minimum. Secondly, fluctuations well beyond the Sachs- Wolfe 
plateau enter the horizon during radiation domination. Unlike adiabatic fluctuations, 
the temperature component of isocurvature fluctuations during this epoch are sup- 
pressed by the radiation density, as can be seen in the denominator of (41). This 
follows since well back into the radiation dominated era the isocurvature condition 
(39) can be satisflcd for flxed magnitude of axion density fluctuation with vanishing 
small temperature fluctuation. This fall off of isocurvature temperature fluctuations 
with multipoles beyond the Sachs- Wolfe plateau is important in considering how well 
such fluctuations can be probed either by a direct measure of the power spectrum or 
search for non-Gaussian components discussed in the next subsection. 

In order to characterize the observable effects of isocurvature ffuctuations it is 
useful to define the ratio of average power in the isocurvature component to the 
average total power in CMBR temperature fiuctuations 



The total root mean square CMBR temperature fiuctuation was measured by COBE 




(48) 
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to be [40] 

{{5T/T)l,)"' ^ 1.1 X 10-^ (49) 

Because of the rapid fall off of isocurvature power, the ratio (48) receives contribu- 
tions mainly from low multipoles on the Sachs- Wolfe plateau. And since, as discussed 
above, axion induced isocurvature temperature fluctuations are fairly well approxi- 
mated by (47) for these low multipoles, the isocurvature power ratio (48) may be 
approximately related to the mean square fluctuations and average value of the ini- 
tial misalignment angle using (43) by 

_ / 6 y (n„/ti,„)' {s'i) . , / 6 \ ^ „ , 20? + 



A full calculation of the relation between a and the variance and average value of the 
initial misalignment angle utilizing the full isocurvature power spectrum illustrated 
in Fig. 1 would not differ signiflcantly from the approximation (50). 

It is instructive to consider the numerical relation (50) between the axion fluctua- 
tion induced isocurvature power fraction and the variance and average initial misalign- 
ment angle in each of the cosmological scenarios for formation of the axion condensate 
discussed in section 2. For the case of smaller values of the Peccei-Quinn scale in which 
the axion condensate forms during a radiation dominated era in the high temperature 
axion potential and without any dilution after its formation, the relic axion density 
is given by (7) in section 2.1. This relic density with C = 0.018 and Aqcd = 200 
MeV and ignoring anharmonic effects, along with the WMAP total matter density 
measurement (45) and COBE total root mean square CMBR temperature fluctuation 
measurement (49) gives a fractional isocurvature power (50) of 

"-'^''"'{ji^T (51) 

For the case of larger Peccei-Quinn scales in which axion condensate forms during a 
radiation dominated era but in effectively the zero temperature axion potential and 
again without any dilution after its formation, the relic axion density is given by 
(9) in section 2.1. This rehc density with /cC(^osc)~^ — 1 and ignoring anharmonic 
effects, along with (45) and (49) gives a fractional isocurvature power (50) of 

Finally, for the case in which the axion condensate forms during a matter dominated 
era with subsequent late decay of the dominating particle with a reheat temperature 
^RH ^ Aqcd, the relic density is given by (16) in section 2.3.1. This rehc density 
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with fc^{Tosc) ^ ^ 1 and ignoring anharmonic effects, along with (45) and (49) gives 
a fractional isocurvature power (50) of 




Note that the fractional isocurvature power is strictly bounded by a < 1. So for a 
string/M-theory axion with fa/N ~ 10^^ GeV the large pre-factors in (51)-(53) imply 
that even just the magnitude of the observed temperature fluctuations (49) places a 
stringent limit on the combination ag{29f + (Jg) of variance and average value of the 
initial misalignment angle in any cosmological scenario for the formation of the axion 
condensate in which the axion exists during inflation. 

The isocurvature power ratio (48) may be measured or bounded by fltting the 
observed CMBR temperature fluctuation power spectrum with a linear combination of 
adiabatic and isocurvature components with some assumptions about the underlying 
cosmological model. Using WMAP data and an analysis in which the tilt of both the 
adiabatic and isocurvature components are allowed to vary over fairly large ranges, 
the constraint a <, 0.4 has been obtained [41]. This analysis is very conservative 
and weakens the bound due to some degeneracies in parameters. The magnitude of 
inflationary induced axion isocurvature fluctuations depends on the Hubble constant 
during inflation. And in most theories of inflation the Hubble constant during inflation 
is by definition a very slowly changing quantity. This results in a nearly scale invariant 
spectrum of axion induced isocurvature temperature fluctuations, niso — 1. If this 
prior were imposed, the resulting bound on a using WMAP data would likely be 
at the many percent level. Since most of the power in isocurvature temperature 
fluctuations is restricted to low multipoles, future measurements will not improve 
knowledge of the power spectrum there since the current measurements are already 
cosmic variance limited even beyond the first peak. However, future measurements 
at higher multipoles along with additional information from polarization, such as 
from PLANCK, will resolve degeneracies in the power spectrum such as the adiabatic 
tilt. And this may well allow the fractional isocurvature power a to be bounded at 
the percent level. It is worth noting that gravity waves also contribute to CMBR 
temperature mainly at low multipoles on the Sachs- Wolfe plateau. If gravity waves 
are indeed detected in a future CMBR experiment, then even less of the budget of 
observed power at low multipoles can be due to a isocurvature component which 
would further strengthen the bound on a. 

3.3 Non-Gaussianity 

Primordial inflationary fluctuations in a QCD axion fleld give rise to CMBR tempera- 
ture fluctuations. These fluctuations have the distinctive feature of being isocurvature 
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[5] with an angular power spectrum which differs from that of adiabatic fluctuations, 
as described in the previous subsection. Axion induced CMBR temperature fluctu- 
ations turn out also to have a non-Gaussian component [6] as reviewed below. This 
additional distinctive feature provides another probe for the existence of a QCD axion 
field during inflation. 

Quantum fluctuations of a free massless fleld generated during inflation are Gaus- 
sian distributed. Slow roll inflation generally requires an inflaton with exceedingly 
small self coupling and so its fluctuations act essentially like a free fleld. Inflaton fluc- 
tuations in most models are therefore Gaussian to a very high degree. The adiabatic 
temperature fluctuations which ultimately result from inflaton induced fluctuations of 
the scalar curvature are then also Gaussian distributed. As discussed above, as long 
as the Peccei-Quinn symmetry remains unbroken, a QCD axion is also essentially a 
free massless field during infiation and so the initial misalignment angle undergoes 
Gaussian distributed fluctuations. However, the axion induced CMBR temperature 
fluctuations are proportional to fluctuations in the axion number density when the 
axion condensate is formed, as discussed in section 3.2. But the number density is 
proportional to the square of the misalignment angle ignoring anharmonic effects, 
Ua (X 9'^, while it is the fluctuations in the misalignment angle, 56, which are Gaus- 
sian distributed. This non-linearity introduces a non-Gaussian component in the 
axion induced temperature fluctuations. 

On large angular scales the axion induced temperature fluctuations (47) may be 
related to fluctuations in the initial misalignment angle through the relation (42) as 



The flrst term in the numerator proportional to 5^ is Gaussian distributed with zero 
mean, while the remaining terms together, [SdY — ag, have a distribution with 
zero mean. So primordial Gaussian distributed fluctuations in the initial misalignment 
angle lead in general to a combination of Gaussian and distributed isocurvature 
temperature fluctuations. In the limit of large average misalignment angle compared 
with the variance, 9i ^ ag, the axion induced temperature fluctuations (54) are 
mostly Gaussian, whereas in the opposite limit, 9i <C ere, they are mostly x^- The 
relative importance of the Gaussian and fluctuations is determined by the ratio 
9i/a0. The temperature fluctuation (54) in these two limits is 



For 9i ^ ag the magnitude of the axion induced temperature fluctuations are sup- 
pressed by the relic axion density and additionally suppressed on average by ag/9i. 
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However, for 6i <^ ag they are only suppressed by the rehc axion density since the 
fractional axion fluctuations are order one in this case. The latter limit turns out to 
be most important in obtaining the most conservative bounds discussed in section 4. 

Gaussian fluctuations with zero mean are completely defined in terms of the two 
point function. All odd point functions vanish, and all even point functions are prod- 
ucts of the two point function. So in order to test for non-Gaussianity it is necessary 
to measure three or higher point functions. In the present context the simplest possi- 
bility to characterize a non-Gaussian component in CMBR temperature fluctuations 
is to measure the average three point function of the temperature fluctuation. This 
is conveniently normalized as a dimensionlcss skcwness 

_ {(ST/Tf) 

[{ST/Tff' 

For a pure Gaussian distribution the skewness vanishes while for a pure distribu- 
tion with zero mean 5*3 = —y/S. Of course the full three point function bi-spectrum 
contains additional information which would be invaluable in discerning the origin 
of non-Gaussianity in CMBR temperature fluctuations if ever observed. However, 
in the absence of such positive measurements the dimensionlcss skewness (56) pro- 
vides a good dimensionlcss normalized measure to characterize the magnitude of 
non-Gaussian components. 

In order to evaluate the dimensionlcss skewness (56) we consider the special case 
in which CMBR temperature fluctuations are an uncorrelated sum of adiabatic plus 
isocurvature contributions 

ST ST ST 

-TfT-^ +^ (57) 

1 -t ad -£ iso 

This is the case if inflaton fluctuations are responsible for the dominant adiabatic 
temperature fluctuations observed in the CMBR, and the only other massless held 
during inflation is a QCD axion. Inclusion of gravity wave tensor B-mode contribu- 
tions from inflationary metric fluctuations would not modify the results below. And 
inclusion of other sources of non-Gaussianity would only allow for a smaller axion 
contribution for a flxed total skewness of CMBR temperature fluctuations, and so 
could only result in stronger bounds than those given in section 4. Since the infla- 
tionary axion fluctuation induced isocurvature temperature fluctuations have both 
Gaussian and components, the total temperature fluctuation (57) may be written 
in terms of a sum of distributions 

ST STi STa ^ ^ 

where (pi parameterizes the Gaussian distributed inflaton induced temperature fluc- 
tuations, (pa the Gaussian distributed component of the axion induced temperature 
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fluctuations, and ipa the distributed non-Gaussian component of the axion in- 
duced temperature fluctuations. Each distribution by definition has zero mean, 
{4>i) = {4>a) = {ipa) = 0. Note that (pa and ipa are not independent since they 
both arise from axion fluctuations whereas 0/ is independent of both 0a and ipa since 
it arises from uncorrelated fluctuations of an independent fleld. So cross correlations 
between 0/ and (pa or ifba vanish. With this, the two and three point temperature 
fluctuation correlation functions from (58) are then 

/5TST\ 1 5Ti 5Ti \ / 5Ta 5Ta 



'6T6T6T\ /6Ta6Ta6Ta 



+ -^-^ = (070/) + {(Pa'Pa) + i^a^a) (59) 



rj^ rj^ rjn / \ rjn rjn rjn 



3 {(paMa) + 3 {(palpalpa) + {ipa'ipa'ipa) (60) 



The three point function of course only depends on the axion fluctuation induced 
isocurvature components since by assumption the adiabatic inflaton fluctuations are 
purely Gaussian distributed. 

Now the dimensionless skewness deflned in (56) is an average over all angular 
scales. Axion isocurvature contributions through the three-point function (60) how- 
ever come mainly from low multipoles on the Sachs- Wolfe plateau, just as for the 
isocurvature power discussed in section 3.2. On these large angular scales the rela- 
tion between the axion induced isocurvature temperature fluctuation and fluctuations 
of the initial misalignment angle is fairly well approximated by (54). The relation 
between the axion induced Gaussian and distributed temperature fluctuation dis- 
tributions (pa and ipa in (58) and the fluctuations in the initial misalignment angle 
may then be approximated by 



6 Qa '^diSe 



Gaussian (61) 
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With this, the three point function (60) and dimensionless skewness (56) for the axion 
induced isocurvature component of temperature fluctuations may be approximately 
related to the mean square fluctuations and average value of the initial misalignment 
angle by 

6\3 {na/^mf A^el + al 



•15^ {{6T/T)iy i(^i+^e) 



A full calculation of the relation between 6*3^180 and the variance and average value of 
the initial misalignment angle utilizing the full isocurvature power spectrum would not 
differ signiflcantly from the approximation (63) . Note that this skewness is necessarily 
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negative as a result of the nature of the non-Gaussian component of axion fluctu- 
ations. This would be the first test that a positive measurement of non-Gaussianity 
in CMBR temperature fiuctuations is in fact due to a QCD axion. 

It is instructive to consider the numerical relation between the axion fiuctuation 
induced isocurvature skewness (63) and the variance and average value of the initial 
misalignment angle in each of the cosmological scenarios for formation of the axion 
condensate discussed in section 2. For the case of smaller values of the Peccei-Quinn 
scale in which the axion condensate forms during a radiation dominated era in the 
high temperature axion potential and without any dilution after its formation, the 
relic axion density is given by (7) in section 2.1. This relic density with C = 0.018 and 
Aqcd = 200 MeV and ignoring anharmonic effects, along with the WMAP total mat- 
ter density measurement (45) and COBE total root mean square CMBR temperature 
fluctuation measurement (49) gives an isocurvature skewness (63) of 



For the case of larger Peccei-Quinn scales in which axion condensate forms during a 
radiation dominated era but in effectively the zero temperature axion potential and 
again without any dilution after its formation, the relic axion density is given by 
(9) in section 2.1. This rehc density with /cC(^osc)~^ — 1 and ignoring anharmonic 
effects, along with (45) and (49) gives an isocurvature skewness (63) of 



Finally, for the case in which the axion condensate forms during a matter dominated 
era with subsequent late decay of the dominating particle with a reheat temperature 
^RH ^ Aqcd, the relic density is given by (16) in section 2.3.1. This relic density 
with fci{Tosc)^^ — 1 and ignoring anharmonic effects, along with (45) and (49) gives 
an isocurvature skewness (63) of 



Note that the axion fiuctuation induced isocurvature skewness is negative definite and 
strictly bounded by ^s^iso > —V8, and as discussed below is already bounded to be 
significantly smaller in magnitude. So for a string/M-theory axion with fa/N ~ 10^^ 
GeV, similar to the case of isocurvature power fraction, the large prefactors in (64)- 
(66) imply that even just the magnitude of the observed temperature fluctuations 
(49) places a stringent limit on the combination (J^(3^f -|- aj) of variance and average 
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value of the initial misalignment angle in any cosmological scenario for the formation 
of the axion condensate in which the axion exists during inflation. 

The dimensionlcss skcwncss (56) may be measured or bounded from the observed 
CMBR temperature fluctuation three point function. As reviewed in appendix C, 
it has become common in the literature to parameterize possible non-Gaussianity 
in terms of a non-linear product of Gaussian-distributed temperature fluctuations 
proportional to a dimensionless parameter /jvl. As derived in (91) of appendix C, 
the dimensionless skewness may be related to this parameter at leading order by 



The WMAP bound on /nl for multipoles £ < 65 is roughly /nl = -150t?g{5 [42]. 
Since the dimensionless skewness for axion induced non-Gaussianity is strictly neg- 
ative, 5'3,iso < 0, we conservatively take Jnl > —300 for this case. With the total 
root mean square CMBR temperature fluctuation (49) measured by COBE [40] this 
bound then corresponds to 5'3^iso ^ — 6 x 10~^. 

Since the magnitude of isocurvature temperature fluctuations falls rapidly be- 
yond the Sachs-Wolfe plateau, as discussed in section 3.2, the inclusion of multi- 
poles much beyond i > 65 would not significantly improve the search for, or bound 
on, a non-Gaussian component in isocurvature temperature fluctuations. And since 
current measurements at these low multipoles arc already cosmic variance limited, 
future measurements will also not signiflcantly improve the current bound on a non- 
Gaussian component in isocurvature temperature fluctuations (although the bound 
on a non-Gaussian component in adiabatic temperature fluctuations, which do not 
fall as rapidly at large multipoles, would be improved). In addition, future reduc- 
tions in degeneracies of parameters describing the dominant adiabatic power spectrum 
from measurements at higher multipoles will also not help improve the bound on a 
non-Gaussian component in isocurvature temperature fluctuations. So unlike the frac- 
tional isocurvature power, a, discussed in section 3.2, the bound on the isocurvature 
skewness, S^^iso is unlikely to improve signiflcantly. 



4 Constraints on the Axion from the CMBR 



It seems unlikely that a string/ M-theory QCD axion with fa/N ~ lO-*^^ GeV will ever 
be probed directly in laboratory experiments. Such an axion can, however, give rise 
to in principle measurable isocurvature [5] and non-Gaussian [6] signatures in CMBR 
temperature fluctuations, as reviewed in section 3. In this section we show that a 
detection of primordial gravity waves interpreted as arising from an early epoch of 
inflation places a lower limit on these observable effects. The current bounds on 
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isocurvature and non-Gaussian components of CMBR temperature fluctuations are 
then shown to be sufficient to rule out a string/M-theory QCD axion if gravity waves 
interpreted as arising from inflation with a Hubble constant i^ing ^ 10^^ GeV are 
observed by future experiments. 

Consider flrst the axion fluctuation induced signature of an isocurvature compo- 
nent in CMBR temperature fluctuations. The goal here is to assess how well CMBR 
experiments can probe a string/M-theory QCD axion independent of any cosmo- 
logical assumptions. This requires considering in turn all the relevant cosmological 
scenarios for formation of the axion condensate presented in section 2. In addition it 
requires considering the smallest possible axion induced isocurvature component in 
CMBR temperature fluctuations consistent with other future CMBR measurements; 
this implies the weakest possible bound. The numerical relations between the frac- 
tional isocurvature power a and variance and average value of the initial misalignment 
angle for the various cosmological scenarios for formation of the axion condensate con- 
sidered in section 2 arc given in (51)-(53) of section 3.2. In each case the fractional 
isocurvature power is proportional to the combination ag{29f + cr^). This is clearly 
minimized in the limit in which the average misalignment angle is much smaller than 
the variance of inflationary fluctuations in the initial misalignment angle, 9"^ ^ aj. In 
this limit of small average misalignment angle, the relic axion density is due mostly 
to inflationary fluctuations of the axion fleld as discussed in section 2.4. And the 
fractional axion fluctuations are order one in this limit, {SD ~ 2 from (43). 

Axion induced isocurvature components in CMBR temperature fluctuations of 
course require primordial fluctuations in the axion fleld. As long as the Peccei-Quinn 
symmetry remains unbroken during inflation, a string/M-theory QCD axion under- 
goes quantum fluctuations during this epoch with magnitude (25) determined by the 
inflationary Hubble constant. As discussed in section 3.1 all massless fields undergo 
quantum fluctuations during inflation including the metric which ultimately gives rise 
to tensor B-modes in CMBR temperature fluctuations. Observation of these modes 
interpreted as arising from inflation would establish the Hubble constant during infla- 
tion and the variance of fluctuations in the axion initial misalignment angle for a given 
Peccei-Quinn scale, as given in (36) in terms of the tensor to scalar ratio. This would 
in turn establish a lower limit for the isocurvature power fraction. So obtaining a 
deflnitive lower limit on an axion induced isocurvature component in CMBR temper- 
ature fluctuations independent of cosmological assumptions requires the observation 
of primordial gravity waves interpreted as arising from inflation. 

The lower limit on a in terms of the Hubble constant during inflation for a given 
cosmological scenario for the formation of the axion condensate may be obtained 
from (51)-(53) in the limit 6f <^ discussed above, along with the relation (36) 
between the variance of the axion initial misalignment angle and inflationary Hubble 
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Figure 2: Minimum allowed axion Peccei-Quinn scale fa/N as a function of iJinf for 
isocurvature temperature fluctuation power fractions a = {{6T /T)f^^) / {{6T /Ty) = 
1, 10~^, 10~^ in the limit that the average initial misahgnment angle is small compared 
to inflation induced fluctuations, 6f <^ (j| = {Hi^i/{2'Kfa/N)Y. For Of > the 
minimum allowed fa/N grows like Of. The dotted lines correspond to no dilution of 
the axion condensate after formation while the solid lines correspond to maximum 
possible dilution by a late decaying particle with Trh — 6 MeV. 

constant. This may also be obtained directly from (50) in the hmit {SD ~ 2 discussed 
above, along with the minimum fluctuation dominated relic axion densities (27)-(29) 
in terms of the inflationary Hubble constant. If gravity wave contributions to CMBR 
temperature fluctuations are in fact observed and interpreted as arising from inflation, 
then the lower limit on a can be inverted to give a minimum allowed Peccei-Quinn 
scale in terms of the implied Hubble constant during inflation and the present bound 
on the isocurvature power fraction. 

The simplest cosmological scenario to consider is obtained if the axion condensate 
is formed during the radiation dominated era without any dilution after its formation. 

In this case for small values of Peccei-Quinn scale the isocurvature power fraction (51) 
with Of <C (tI and the relation (36) yields a lower limit on the Peccei-Quinn scale as 
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outlined above in terms of a and the Hubble constant during inflation of 



/./^ ^ ^ «eV e^<al (68) 

In the same case for large values of Peccei-Quinn scale, the isocurvature power fraction 
(52) with 9l -C al and the relation (36) yields a lower hmit on the Peccei-Quinn scale 

These bounds on the Peccei-Quinn scale are plotted in Fig. 2 as the diagonal dotted 
and sohd lines for a — 1, 10~^, 10~^. The bounds for large and small Peccei-Quinn 
scale coincide for fa/N ~ 6 x 10^^ GeV. As discussed in section 2.1 for Peccei-Quinn 

scales in this transition region the relic axion density calculation suffers unknown 
strong QCD uncertainties but should asymptotically approach the limiting expres- 
sions for scales well outside the transition region. 

In the remaining cosmological scenario to consider, the axion condensate is formed 
during a matter dominated era with subsequent decay of the dominating particle. In 
this case with <^ a^, the maximum possible dilution consistent with BBN discussed 
in section 2.3.1, Trh — 6 MeV, and the relation (36) between the variance of the 
initial misalignment angle and Hubble constant during inflation, the isocurvature 
power fraction (53) is bounded by 

Note that this bound is independent of the Peccei-Quinn scale since the explicit 
(fa/N)^ dependence in (53) cancels the {fa/N)~^ dependence in a^. This can also 
be seen directly in (50) with {So) ~ 2 as discussed above in this limit, along with 
the minimum inflationary fluctuation induced relic axion density (29) which is inde- 
pendent of Peccei-Quinn scale as discussed in section 2.4. This implies that if the 
bound (70) is not satisfied, then for an axion which exists during inflation, all Peccei- 
Quinn scales arc cxchidcd for which the axion condensate forms during the matter 
dominated era before subsequent decay of the dominating particle. As discussed in 
section 2.4, this occurs with maximum dilution, Trh ~ 6 MeV, for Peccei-Quinn 
scales fa/N < 2 x 10^° GeV. For larger Peccei-Quinn scales the axion condensate 
forms during the subsequent radiation dominated era after reheating from decay of the 
dominating particle. So for fa/N > 2 x 10^° GeV the bound (69) is recovered in any 
scenario for formation of the axion condensate. The bounds obtained from (70) for 
a = 1, 10^^, 10^^ with maximum dilution, Trh ~ 6 MeV, are plotted as vertical solid 
lines in Fig. 2. Since the axion condensate forms in essentially the zero temperature 
potential in this case, these bounds do not suffer large QCD uncertainties. 
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Now consider the axion fluctuation induced signature of a non-Gaussian compo- 
nent in CMBR temperature fluctuations. The numerical relations between the axion 
induced isocurvature skewness, S^^iso, and variance and average value of the initial 
misalignment angle for the various cosmological scenarios for formation of the ax- 
ion condensate considered in section 2 are given in (64)-(66) of section 3.3. In each 
case the skewness is proportional to the combination ag{36f + ag). Just as for the 
isocurvature power fraction, this combination is minimized for 6"^ <^ aj which yields 
in turn the weakest possible bound. In this limit of small initial misalignment angle, 
the relic axion density is due mostly to inflationary fluctuations of the axion fleld as 
discussed in section 2.4. And the axion induced isocurvature temperature fluctua- 
tions are approximately distributed as discussed in section 3.3. The lower limit 
on 5*3^180 may be obtained in the 6f <^ (Tq limit from (64)-(66) along with the relation 
(36) between the variance of the axion initial misalignment angle and inflationary 
Hubble scale. This may also be obtained directly from (63) in this limit along with 
the minimum fluctuation dominated relic axion densities (27)-(29) in terms of the 
inflationary Hubble constant. Just as for the isocurvature power fraction, if gravity 
wave contributions to CMBR temperature fluctuations are in fact observed and in- 
terpreted as arising from inflation, then the lower limit on 5*3^180 can be inverted to 
give a minimum allowed Peccei-Quinn scale in terms of the implied Hubble constant 
during inflation and the present bound on the isocurvature skewness. 

The simplest cosmological scenario for formation of the axion condensate is that it 
formed during the radiation dominated era without any dilution after its formation. In 
this case for small values of Peccei-Quinn scale the isocurvature skewness (64) yields 
a lower limit of the Peccei-Quinn scale in terms of S^^iso and the Hubble constant 
during inflation, as described above, of 

In the same case for large values of Peccei-Quinn scale, the isocurvature skewness 
(65) with Of <^ (Tq and the relation (36) yields a lower limit on the Peccei-Quinn scale 
of 

These bounds on the Peccei-Quinn scale arc plotted in Fig. 3 as the diagonal dotted 
and solid lines for —S^^iso = 10"^, 10^^. The bounds for large and small Peccei-Quinn 
scale coincide for fa/N ~ 6 x 10^'' GeV, and suffer unknown QCD uncertainties in 
the transition region. 

In the remaining cosmological scenario to consider, the axion condensate is formed 
during a matter dominated era with subsequent decay of the dominating particle. In 
this case with 9'f <^ cr^, the maximum possible dilution consistent with BBN discussed 
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Figure 3: Minimum allowed axion Peccei-Quinn scale /a/iV as a function of H-^ai 
for isocurvature temperature fluctuation skewness S'a = {{5T /T)l^^ / {{5T /TYY^'^ — 
— IQ-i, — 10~^ in the limit that the average initial misalignment angle is small com- 
pared with inflation induced fluctuations, Of — {Hi^i/{2'Kfa/N)Y. For 9f > 
the minimum allowed fa/ N grows like Of^^. The dotted lines correspond to no dilution 
of the axion condensate after formation while the solid lines correspond to maximum 
possible dilution by a late decaying particle with Trh — 6 MeV. 



in section 2.3.1, Trh — 6 MeV, and the relation (36) between the variance of the initial 
misalignment angle and Hubble constant during inflation, the isocurvature skewness 
(66) is bounded by 

Just as for the isocurvature power fraction in this scenario for formation of the axion 
condensate, this bound is independent of the Peccei-Quinn scale since here the explicit 
{fa/NY dependence in (66) cancels the (fa/N)'^ dependence in a^. This imphes 
that if the bound (73) is not satisfied, then for an axion which exists during inflation, 
all Pcccci-Quinn scales arc cxchidcd for which the axion condensate forms during 
the matter dominated era before subsequent decay of the dominating particle. As 
discussed in section 2.4, this occurs with maximum dilution, Trh ~ 6 MeV, for 
Peccei-Quinn scales fa/N < 2 x 10^° GeV. For larger Peccei-Quinn scales the axion 
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condensate forms during the subsequent radiation dominated era after reheating from 
decay of the dominating particle. So for fa/N > 2 x 10^*^ GeV the bound (72) is 
recovered in any scenario for formation of the axion condensate. The bounds obtained 
from (73) for —5*3,180 = 10~^, 10~^ with maximum dilution, Trh ^ 6 MeV, are plotted 
as vertical solid lines in Fig. 3. Since the axion condensate forms in essentially 
the zero temperature potential in this case, these bounds do not suffer large QCD 
uncertainties. 

The bounds arising from isocurvature and non-Gaussian contributions to CMBR 
temperature fluctuations presented in Figs. 2 and 3 on the Peccei-Quinn scale for an 
axion which exists during inflation are quite severe if primordial gravity waves are 
observed in future experiments. As discussed in section 3.1, the search for tensor B- 
mode components in CMBR temperature fluctuations with the PLANCK polarimetry 
experiment may reach an ultimate sensitivity to the Hubble constant during inflation 
of -ffinf > 5 X 10^^ GcV [33, 34], with a likely benchmark for a positive observa- 
tion of i^inf > 10^"^ GeV. Even with the rather conservative current bound on the 
isocurvature power fraction discussed in section 3.2 of a > 0.4 [41], if the PLANCK 
experiment did in fact observe tensor B-modes interpreted as arising from inflation 
with iJinf > 10^'^ GeV, then (69) would require that the Peccei-Quinn scale be 
bounded by at least fa/N ^ 6 x 10^^ GeV. And in the same case the current bound 
on the isocurvature skewness discussed in section 3.3 of —S^ ^ 6 x 10^^ [42] with (72) 
would require that the Peccei-Quinn scale be bounded by at least fa/N > 3 x 10^° 
GeV. If CMBR temperature fluctuation tensor B-modes were discovered with mag- 
nitude not too far below the current WMAP implied limit on the Hubble constant 
during inflation discussed in section 3.1 of /Jjnf > 2.8 x 10"*^^ GeV [15] the bounds 
would be even more severe. For example, if it turned out that ifinf = 2 x 10^^ GeV the 
current isocurvature power fraction bound would require fa/N > 10^^ GeV, while 
the current isocurvature skewness bound would require fa/N > 5 x 10^^ GeV. 

If primordial gravity waves are in fact observed, the bound on the Peccei-Quinn 
scale from the present bound on a non-Gaussian component in CMBR temperature 
fluctuations would be somewhat more stringent than that from the present bound on 
the isocurvature power fraction. However, as discussed at the end of sections 3.2 and 
3.3, the bound on isocurvature non-Gaussianity is not expected to improve in the 
future, while the bound on isocurvature power fraction should improve considerably 
either through a re-analysis of current WMAP data and/or future CMBR temperature 
fluctuation data at larger multipoles. So ultimately the bound from the isocurvature 
power fraction would likely be somewhat better. In comparing the isocurvature power 
fraction with isocurvature skewness bounds note that the parametric dependence on 
the Hubble constant during inflation of (68) and (71) for low Peccei-Quinn scales 
are identical, as are those of (69) and (72) for high Peccei-Quinn scales. In the case 
that the axion condensate is formed during a matter dominated era, the parametric 
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dependence of the bounds (70) and (73) on the Hubble constant during inflation and 
reheat temperature after the axion condensate forms are different however. 

For lower values of the Hubble constant during inflation the implied bounds on the 
Peccei-Quinn scale given in Figs. 2 and 3 are reduced. But as long as i!finf ^ 10^'^ 
GeV, even with the current bounds on the CMBR isocurvature power fraction and 
isocurvature skewness, the bounds imply that in this case the only allowed scenario 
for formation of the axion condensate would be during the radiation dominated era 
at a temperature lower than minimum allowed reheat temperature of Trh ~ 6 MeV 
implied by consistency with BBN. In this regime the axion condensate forms in essen- 
tially the zero temperature potential and the relic density calculation does not suffer 
any strong QCD uncertainties. So the bounds presented here that would apply with a 
positive measurement of primordial gravity waves by PLANCK are largely free from 
such uncertainties. 

The bounds discussed above on the Peccei-Quinn scale which would result if pri- 
mordial gravity waves are observed in a future experiment arc quite severe. However, 
it should be noted again that these bounds arc the weakest possible given all scenarios 
for formation of the axion condensate, and average initial misalignment angle. These 
weakest possible bounds are obtained in only an extraordinarily tiny sliver of the 
full parameter space of initial misalignment angles. Outside these tiny slivers even 
stronger bounds would apply. For example, even for the minimum Hubble constant 
during inflation of roughly ifinf ~ 10^'^ GcV to which the PLANCK polarimetry ex- 
periment will be sensitive for a positive measurement, the bound on the Peccei-Quinn 
scale quoted above of fa/N > 3 x lO^'^ GeV which would result from the current 
bound on the isocurvature skewness, corresponds to an initial misalignment angle 
of at most \9i\ < 5 x 10~^. For fixed i?inf and isocurvature skewness, the bound 
on the Peccei-Quinn scale grows like Of^^ for 6f > a^, and for fixed isocurvature 
power fraction like 9f. This also implies that even if the bounds were just barely 
saturated, relic axions would comprise a very small fraction of the CDM. This can 
be seen directly, for example, in the relation (50) between the isocurvature power 
fraction, fractional axion fiuctuations, relic axion density, and magnitude of the total 
mean square temperature fluctuations. Since as discussed above, {SD ~ 2 results in 
the weakest possible bound, the current bound of a < 0.4 in this limit, along with 
the COBE measurement (49) of the magnitude of temperature fluctuations, implies 
that if the bounds are just saturated then ^ 1 x 10~^ at least. 

Finally, let us consider what implications the stringent bounds on the Peccei- 
Quinn scale given above, if realized by a future measurement of primordial gravity 
waves, would have on a string/M-theory QCD axion. First, since a string/ M-theory 
QCD axion exists during inflation these bounds would apply. In all cases for an 
observation of gravity waves which implied a Hubble constant during inflation of 
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Hini ^ 10^^ GeV, definitive bounds on the Peccei-Quinn scale would result as dis- 
cussed above. For this Hubble constant during inflation, the current bound on isocur- 
vature power fraction would require at least fa/N > 6 x 10^^ GeV and the current 
bound on the isocurvature skewness would require at least fa/N > 3 x 10^° GeV. 
These bounds would exceed the four-dimensional Planck scale of Mp ~ 2.4 x 10^^ GeV. 
Now obtaining a Peccei-Quinn scale which is in excess of the four-dimensional Planck 
scale seems problematic in a fundamental theory of gravity. The maximum scale for 
compact moduli in string/M-theory is strongly believed to be the four-dimensional 
Planck scale [14]; no counter example is known. With this restriction we are led to the 
strong conclusion that if primordial gravity waves interpreted as arising from inflation 
are observed with an implied Hubble constant iJjnf > 10^^ GeV, which includes the 
range of expected sensitivity of the PLANCK CMBR polarimetry experiment, then a 
string/M-theory compact modulus can not play the role of the QCD axion and solve 
the strong CP problem. 



5 Remaining Axion Windows 

The bounds presented in section 4 rather convincingly require a QCD axion which 
exists during inflation to have a Peccei-Quinn scale of at least fa/N ^ 3 x 10^*^ GeV if 
gravity waves interpreted as arising from inflation with Hubble constant i^inf ^ 10^^ 
GeV are observed by the PLANCK polarimetry experiment, and at least fa/N ^ 
5 X 10^^ GeV if gravity waves are observed just below the current WMAP bound. 
As discussed above, these bounds are the weakest possible considering all relevant 
cosmological scenarios for formation of the axion condensate and all possible average 
initial misalignment angles, and are therefore quite conservative. However, it is worth 
discussing under what conditions these bounds do not apply and possible loopholes 
for evading the bounds. 

In some scenarios for realizing the Peccei-Quinn mechanism, the QCD axion arises 
from spontaneous breaking of an anomalous global U{1) symmetry. In this case 
the Peccei-Quinn phase transition to the broken phase could in principle take place 
after inflation. The axion fleld would then not exist during inflation, so the CMBR 
temperature fluctuations bounds discussed here, which depend on the existence of 
the axion during inflation, would not apply. However, in such a scenario global 
axion cosmic strings are formed during the Peccei-Quinn phase transition by the 
Kibble mechanism. These cosmic strings reach a scaling distribution through axion 
radiation and chopping off string loops. The cosmic string radiated axions ultimately 
contribute to the relic axion density. The total number density of string radiated 
relic axions can conservatively be estimated to be roughly an order of magnitude 
more than those arising from coherent production in the zero mode [44] . And in this 
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scenario since the QCD vacuum angle winds around the axion strings 9 G (0,27rA^], 
the average initial misalignment angle squared is {9"^) = 7r^/3. So without a dilution 
from a late decaying particle, the bound (11) discussed in section 2.1, including cosmic 
string radiated axions, would be strengthened to roughly fa/N < few x 10^*^ GeV. 
With a late decaying particle with the lowest possible reheat temperature consistent 
with BBN constraints, Trh 6 McV, the bound (17) discussed in section 2.3 would 
likewise be strengthened to at best fa/N < 10^'' GeV. The lower bound on the 
Peccei-Quinn scale in any scenario for realizing the Peccei-Quinn mechanism with a 
QCD axion comes from the coohng rate of supernova SN1987A which requires roughly 
fa/N <; 10^° GeV [45]. These upper cosmological bounds and lower astrophysical 
bound define the open window for a QCD axion in a cosmological scenario with a 
Peccei-Quinn phase transition after inflation. 

In scenarios in which the QCD axion does indeed exist during inflation, the most 
important consideration for applicability of the bounds discussed in section 4 is the 
observation of gravity waves either from tenor B-mode contributions to CMBR tem- 
perature fluctuations or other means. If these are not observed then the allowed 
axion window can be quite wide with in principle no cosmological upper bound. In 
this case the Hubble constant during inflation might in fact be unobservably small. 
Searches for isocurvature or non-Gaussian components of CMBR temperature fluc- 
tuations then do not provide definitive probes for a QCD axion which exists during 
inflation. However, without dilution of the axion condensate after formation, exceed- 
ing the bound (11) of fa/N < 3 x 10^^ GeV does requires a small misahgnment angle. 
And as discussed above this is not possible in a scenario with a Pcccci-Quinn phase 
transition after inflation, and so would require that the axion exist during inflation 
and happen to have a small initial misalignment angle. Even with a late decaying 
particle with the lowest possible reheat temperature consistent with BBN constraints, 
Trh ~ 6 MeV, exceeding the bound (17) of fa/N < 3 x 10^^ GeV also requires a 
small misalignment angle, and so could only be realized if the axion existed during 
inflation. 

It is also worth considering whether there are any possible loop holes to the bounds 
presented in section 4 in the case that gravity waves interpreted as arising from in- 
flation are in fact observed in future experiments. As emphasized, application of the 
CMBR temperature fluctuation isocurvature and/or non-Gaussian bounds requires 
the existence of a QCD axion during inflation which undergoes quantum fluctuations, 
which in turn requires that the Peccei-Quinn symmetry remain unbroken during in- 
flation. It is possible in principle, however, that the Peccei-Quinn symmetry could 
be broken during inflation, but restored after inflation. In this case the axion would 
possess a non-trivial potential during inflation which could in principle suppress quan- 
tum fluctuations if 'ma,mi ^ -f^inf , where ma,inf is the mass scale of the induced axion 
potential during inflation. Such a scenario requires a signiflcant breaking of Peccei- 
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Quinn symmetry during inflation. This can be obtained in field theory models for 
the origin of a QCD axion designed specifically for this purpose [46]. But for large 
Peccei- Quinn scales the minimum of the axion potential during infiation would have 
to be quite close to the zero temperature minimum in order for the relic axion density 
to be consistent with the observed CDM density. 

For a string/M-theory QCD axion it is not clear how to obtain such a large 
violation of Peccei-Quinn symmetry in the early universe, consistent with the req- 
uisite exceedingly small violation in the current epoch for a successful solution of 
the strong CP problem. One possibility might be a spontaneous breaking of elec- 
troweak SU{2)l x U{1)y gauge symmetry during infiation by a large expectation 
value along the HuH^ D-Hai direction. If the expectation value were sufficiently large 
that all the resulting standard model quark masses where larger than the Hubble 
constant, niq > i^inf, then an appropriately small dilaton expectation value during 
infiation could in principle lead to gluino condensation for SU{3)c during infiation 
with AQCD,inf 3> -ffinf- If in addition symmetry were maximally broken during 

infiation then a gluino mass of at most ms < ifjnf could be generated [43] . If all these 
conditions were met, then mixing of the string/M-theory QCD axion with the phase 
of the gluino condensate would lead to a potential for the axion during infiation with 
mass scale m^ j^j ~ ■'^sAqcD.mf/l/a/^)^- For AqcD.inf sufficiently large, 
might be possible which would suppress quantum fiuctuations. This seems the only 
reasonable scenario for suppressing quantum fiuctuations of a string/M-theory QCD 
axion during infiation. While rather contrived, it may be kept in mind as a possible 
loop hole to the bounds discussed here in the case that primordial gravity waves are in 
fact observed, and considered in more detail at that time. However, it does suffer the 
problem that during infiation the sum of the phases of the Higgs condensate HuH^ 
and the phase of the U{l)ji breaking gluino mass 777.3 would have to be quite close to 
the sum of these phases at zero temperature. Otherwise the minimum of the axion 
potential during infiation would not be close to the zero temperature minimum, and 
would lead to a relic axion density in excess of the observed CDM density. 

Another possible loophole is the form of a late entropy release which could dilute 
the axion condensate after formation. The general case of dilution by a late decaying 
particle was considered in section 2.3.1, and the largest such dilution allowed by the 
successful predictions of BBN was included in the bounds presented in section 4. An 
arbitrary entropy release at a temperature below that at which the axion condensate 
forms, Tosc ^ Aqcd, but with reheat temperature well in excess of Trh <^ 1 MeV 
could leave BBN unmodified while diluting the axion condensate by an arbitrary 
amount. This might in principle occur from a strongly first order phase transition in 
this temperature range. However, such a strong phase transition would require new 
physics at this energy scale with significant tree-level coupling to standard model 
fields in order to achieve successful reheating. This is certainly ruled out by the non- 
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observation of such physics in laboratory experiments. An extreme version of such 
an entropy release would be a very late inflation with Hubble constant H > rUa, 
but with a reheat temperature Trh ^ 1 MeV [47]. A reheat temperature this large 
with such a low Hubble constant would again require significant tree-level coupling 
of the inflaton sector at the minimum of its potential to standard model fields at this 
energy scale, and is also certainly ruled out by laboratory experiments. Even if not 
already ruled out, such scenarios would have to cope with the possible problem of 
over dilution or production of the baryon asymmetry. So an arbitrary entropy release 
after formation of the axion condensate does not seem to be a credible loop hole. 



6 Conclusion 

String/M-theory vacua generically have compact moduli which can potentially act as 
the QCD axion and provide an elegant solution of the strong CP problem. In large 
classes of vacua, such as with unification of the gauge couplings by four- dimensional 
renormalization group running, the Peccei-Quinn scale for a string/ M-theory QCD 
axion can be large, fa/N ~ 10^^ GeV. Obtaining an acceptable cosmological scenario 
to accommodate such an axions can be problematic. In order for relic axion produc- 
tion to be consistent with measurement of the present CDM density, a small initial 
misalignment angle to minimize the relic axion density, and possibly an late entropy 
release to dilute it further, are required. In addition, the late decays of relic saxions 
and axinos present problems in many scenarios for supersymmetry breaking, which 
must be addressed. But it is certainly possible to construct cosmological scenarios 
with a string/M-theory QCD axion in which all the requisite conditions are met and 
problems avoided. 

The goal here was to assess how well precision cosmological measurements could 
probe a string/M-theory QCD axion in any cosmological scenario which can accom- 
modate such an axion. This required considering all possible scenarios which could 
in principle minimize the observational effects of relic axions, consistent with possi- 
ble future precision cosmological measurements. Establishing a minimal magnitude 
for observational effects requires an independent measurement which would imply a 
lower limit on the density of relic axions. Such a measurement would be provided by 
the observation of primordial gravity waves interpreted as arising from infiation. A 
measurement of the power in primordial gravity waves would establish the magnitude 
of inflationary induced quantum fluctuations of the axion fleld. For a given Peccei- 
Quinn scale this would in turn establish a minimum density in the axion condensate 
when formed, arising from at least the inflationary axion fluctuations. Minimizing the 
relic axion density resulting from this formation requires in turn considering a scenario 
with an entropy release by a late decaying particle which dilutes the axion condensate 
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by the maximally allowed amount consistent with the successful predictions of BBN. 

A measurement of the power in primordial gravity waves would also establish 
the minimal contribution of relic axion fluctuations to both isocurvaturc and non- 
Gaussian components in CMBR temperature fluctuations. Bounds already exist from 
WMAP on iso curvature and non- Gaussian contributions to CMBR temperature fluc- 
tuations. Given even these current bounds, we have shown that an observation of 
primordial gravity waves interpreted as arising from an early epoch of inflation with 
Hubble constant iJi^f ^ 10^^ GeV would imply for a QCD axion which exists dur- 
ing inflation, a bound on the Peccei-Quinn scale which is in substantial excess of 
the four-dimensional Planck scale. However, the Peccei-Quinn scale for a string/M- 
theory compact modulus is strongly believed to be bounded from above by the four- 
dimensional Planck scale [14]. The above range of power in primordial gravity waves 
is expected to be within the reach of sensitivity of the PLANCK polarimetry exper- 
iment to tensor B-mode contributions to CMBR temperature fluctuations. In the 
future, space based gravity wave interferometers may also reach this level of sensitiv- 
ity in direct searches for background primordial gravity waves. All this leads to the 
following strong conclusion: 

A string/M-theory compact modulus can not play the role of a QCD axion 
and solve the strong CP problem if gravitational waves interpreted as aris- 
ing from inflation with a Hubble constant ifinf > 10-^^ GeV are observed 
in future experiments. 

The sensitivity to isocurvaturc components in CMBR temperature fluctuations are 
expected to improve considerably beyond the current bound. This implies the possi- 
bility for an interesting corollary to the above conclusion: 

// isocurvaturc and/or non-Gaussian contributions to CMBR tempera- 
ture fluctuations are observed and interpreted as coming from a string/M- 
theory QCD axion, then primordial gravity waves arising from an early 
epoch of inflation with ifjnf ^ 10^^ GeV will not be observed. 

These conclusions require only the assumption that the Peccei-Quinn symmetry re- 
mains unbroken in the early universe. With this mild assumption the bounds which 
imply these conclusions are quite stringent even though they are the weakest possible 
considering all relevant scenarios for the formation of the axion condensate and initial 
misalignment angle. 

The stringent bounds on the Peccei-Quinn scale which would result from a positive 

measurement of primordial gravity waves would apply generally to any QCD axion 
model in which the axion exists during inflation. And the bounds would imply that 
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any such axion with a Peccei-Quinn scale which exceeded the imphed bounds would 
comprise at most a very small fraction of the CDM. 

The coming generation of precision cosmological observations will provide a wealth 
of data on cosmological evolution. This data can have an direct impact on our under- 
standing of very early universe cosmology. As demonstrated here, these observations 
can also have a direct impact on our understanding of fundamental physics at the 
highest possible energy scales. 
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A The Model Independent String Axion 

All perturbative string vacua posses a model-independent axion supermultiplet which 
couples universally to all gauge kinetic terms. The relation between the Peccei- 
Quinn scale, four-dimensional Planck scale, and unified gauge coupling for this axion 
may be determined from the form of the Lagrangian. The gauge kinetic terms after 
compactification to four dimensions and field redefinitions to supersymmetric Einstein 
frame are of the form 



/ 



d^e ]s w^wi + h.c. D + i^fi^^f"'"' (74) 



where is the superfield strength, F^^ the field strength, F^^ — ^ef^^paF"-'"^ the 
dual field strength, and 

S=\ + ^, (75) 

the bosonic component of the axion supermultiplet. The real part of the axion mul- 
tiplet is related to the unified value of the gauge coupling by 
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where au — Qu/^'^- And with the normahzation (74) the ^-angle, 9 — 87r^Im S is 
periodic mod 27r. 

Starting from any of the ten-dimensional string theories, the axion multiplct La- 
grangian kinetic terms in supersymmetric Einstein frame after compactification to 
four dimensions are [48] 



where Mp is the four-dimensional reduced Planck mass, Mp = rrip/vSn. In any model 
the vacuum angle is related to the canonically normalized axion field and Peccei-Quinn 
scale by 

ft- ° 

where is the axion anomaly coefficient. The axion Lagrangian kinetic term is then 

ld,ad^a^^^d,9d^9 (77) 

For the model independent string axion = 1 for gauge groups realized at the first 
Kac-Moody level. Comparing (76) and (77) the Pcccci-Quinn scale is then related to 
the four-dimensional Planck mass and unified gauge coupling by 

/„ = ^/2 ^Mp ~ 10^^ GeV (78) 

4:71 

Previous discussions [4, 49, 26] of the model independent axion Peccei-Quinn scale 
differ from this definition by a factor proportional to g^'^ because of improperly nor- 
malized ^-angle. 



B Zero and Finite Temperature Axion Mass 

At zero temperature the physical axion mass results mainly from the small pion 
component induced by mixing with the QCD pseudo-Goldstone multiplet. This may 
be calculated by either current algebra [50] or chiral Lagrangian techniques with the 
result 

"^«"TT^7^"'-'^'° ^^^^ ^''^ 

where z = niu/md ^ 0.56 is the ratio of up to down current quark masses, and /^r — 93 
MeV is the pion decay constant. The leading corrections to (79) are 0(m^/m|^) from 
Kaon mixing. 
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At finite temperatures above the chiral symmetry breaking scale, T > Aqcd, 
the axion mass may be determined from the ^qco-dependence of the thermal free 
energy. At sufficiently high temperature, T ^ Aqcd, this may be reliably related 
in the dilute instanton gas approximation to the single instanton amplitude. In this 
approximation the ^qco-dependent free energy density resulting from summing over 
both instantons and anti-instantons is related to the integral over instanton density 
n{p) at size p by [20]* 

F(e, T) ~ -2 cos(^) j dp n{p, T) (80) 

For SU (Nc) gauge group with Nf light quark fiavors the finite temperature instanton 
density in this approximation is [20] 



Cn^ f 



2Nc 



n{p, T)^^ j (Xp)^^(det m) e-'^'/^'^'^^e-f^'^^^ (81) 

where Cat^ = (0.260 156)x"^^^"^V[(^c - l)!(^c - 2)!], x = 1-338 76, m is the light 
quark mass matrix, and 



fix) = y (2Ar, + Nf)x'' + 12A{x) ^ 
where A{x) is very well fit by the expression 



1 + -(N,-Nf) 



(82) 



A{x) ~ In 1^1 + yx'j +a[l+ 7(7rx)-^/'] * (83) 

where a — 0.012 897 64 and 7 = 0.158 58. The next to leading order two-loop 
renormalization group improvement of the scale dependent gauge coupling, g'^{p), 
may be related to Aqcd, defined as the pole of the inverse gauge coupling defined at 
this order, by [20] 

^""^ Mn(l/pAQCD) + ^lnln(l/pAQCD) (84) 



where h = |(llA^'c ^ 2A^/) is the coefficient of the one-loop gauge /3-function, and 
62 = |[17A^^ — A^/(13A^^ — 3)/2A^c]- At next to leading order in the instanton density 
(81) the two- loop improvement for the scale dependence of the gauge coupling (84) 
must be included in the instanton amplitude e^^'^ 1^ '^''\ but the scale dependence of 
the At^"^ / g'^{p) prefactor need only be evaluated at the one- loop level. 

*Expression (5.6) of Ref. [20] for the ^-dependence of the free energy should contain a factor of 
2 for summing over both instantons and anti-instantons. 
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With the normahzations given above, the finite temperature axion mass squared, 
F(a,T) -ml(T){fjNfcos[a/{fjN)] = ^ml{T)a' + ■■■ (85) 
in the dilute instanton gas approximation at next to leading order is then 

by""' KcB detm /Aqcd\^^+^-^ 

2Nc-b2/b 

e-^(^) (86) 

where x — pT. In the high temperature limit the integral is dominated by instantons 
of size roughly of order p ~ vr/T. The integral in (86) has been evaluated numerically 
for the physically relevant case of A'c = 3 and various Nf with the result that the 
temperature dependent axion mass may be fit by the functional form [17] 

For the normalizations given above and correctly including both the instanton and 
anti-instanton contributions with Nf = 3 light quark flavors relevant to the oscilla- 
tions temperatures which result for fa/N ~ 10^^ GeV, the numerical parameters of 
the fit (87) are C ~ 0.018, a = |, ^9 ~ 4.0, and d ~ 1.2 [17]. For oscillation tem- 
peratures not too much larger than Aqcd the ln(AQCD/7") term coming from next to 
leading order effects may be neglected. This approximation is used in deriving the 
oscillation temperature given in section 2. 
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C Parameterizing Non-Gaussian Fluctuations 

Scalar temperature fluctuations in the CMBR may generally be related to a gauge 
invariant potential function by 6T/T = $/3. One manner in which non-Gaussian 
temperature fluctuations could arise is that the potential $ is a sum of terms which 
are linear and quadratic in an underlying Gaussian distributed potential function. As 
discussed in section 3.3 temperature fluctuations induced by axion fluctuations are 
in fact of this form. It has become common in the literature to parameterize this 
possibility by a sum of linear and non-linear contributions to the total potential 

$ = $L + /ArL($i-($i)) (88) 

where $l is the underlying Gaussian distributed potential function, and Jnl is a 
constant parameterizing the strength of the non- linear fluctuations quadratic in 
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This parameterization of non-Gaussian fluctuations may be related to the dimen- 
sionless skewness (56) defined in section 3.3. To see this, first consider the mean 
square variance of the temperature fluctuation with the parameterization 



{^T/T f) = ^ (l + 2/^,4 J - ^ (89) 

where — ($|) and the second approximate equal sign results if the non-Gaussian 

non-linear contributions to the potential are small compared with the linear Gaussian 
contributions, Jml'^^l ^ ^' ^^^^ consider the mean cubic variance of the temper- 
ature fluctuations with the parameterization (88). To leading order in f^i it is a 
product of two linear Gaussian contributions and one non-linear quadratic Gaussian 
contribution to the potential [c.f. the leading term on the right hand side of (60)] 

{ST/Tf) I ($i - = ^/iVL< (90) 

where ($^) = 3(t|,^. The dimensionless skewness (56) is then related to Jnl at leading 
order by 

= ^^^^'^]l/2 ^ ^fNLa^, ^ l^fNL ({ST/T)iy (91) 



{6T/T) 



References 

[1] S. Dimopoulos and H. Gcorgi, "Softly Broken Supcrsymmetry And SU (5)," Nucl. 
Phys. B 193, 150 (1981); S. Dimopoulos, S. Raby and F. Wilczek, "Supcrsym- 
metry And The Scale Of Unification," Phys. Rev. D 24, 1681 (1981). 

[2] N. Kalopcr, M. Kleban, A. E. Lawrence and S. Shenker, "Signatures of short 
distance physics in the cosmic microwave background," Phys. Rev. D 66, 123510 
(2002), hep-th/0201158. 

[3] R. D. Peccei and H. R. Quinn, "CP Conservation In The Presence Of Instantons," 
Phys. Rev. Lett. 38, 1440 (1977); S. Weinberg, "A New Light Boson?," Phys. 
Rev. Lett. 40, 223 (1978); F. Wilczek, "Problem Of Strong P And T Invariance 
In The Presence Of Instantons," Phys. Rev. Lett. 40, 279 (1978). 

[4] T. Banks and M. Dine, "The cosmology of string theoretic axions," Nucl. Phys. 
B 505, 445 (1997), hep-th/9608197. 

[5] D. Seckel and M. S. Turner, "'Isothermal' Density Perturbations In An Axion 
Dominated Inflationary Universe," Phys. Rev. D 32, 3178 (1985); A. D. Linde, 



47 



"Generation Of Isothermal Density Perturbations In The Inflationary Universe," 
Phys. Lett. B 158, 375 (1985). 

[6] D. H. Lyth and E. D. Stewart, "Constraining the inflationary energy scale from 
axion cosmology," Phys. Lett. B 283, 189 (1992); D. H. Lyth, "Axions and 
inflation: Vacuum fluctuations," Phys. Rev. D 45, 3394 (1992). 

[7] http: / / astro.estec.esa.nl/SA-general/Projects/Planck/ 

[8] D. H. Lyth, "A Limit On The Inflationary Energy Density From Axion Isocur- 
vature Fluctuations," Phys. Lett. B 236, 408 (1990). 

[9] A. D. Linde, "Axions in inflationary cosmology," Phys. Lett. B 259, 38 (1991). 

[10] M. S. Turner and F. Wilczek, "Inflationary Axion Cosmology," Phys. Rev. Lett. 
66, 5 (1991). 

[11] E. P. S. Shellard and R. A. Battye, "Cosmic axions," astro-ph/9802216; 
E. P. S. Shellard and R. A. Battye, "On the origin of dark matter axions," 
Phys. Rept. 307, 227 (1998), astro-ph/9808220. 

[12] S. D. Burns, "Isentropic and isocurvature axion perturbations in inflationary 
cosmology," astro-ph/9711303. 

[13] M. Kawasaki and T. Yanagida, "Are isocurvature fluctuations of the M-theory 
axion observable?," Prog. Theor. Phys. 97, 809 (1997) hep-ph/9703261. 

[14] T. Banks, M. Dine, P. J. Fox and E. Gorbatov, "On the possibility of large axion 
decay constants," JCAP 0306, 001 (2003), hep-th/0303252. 

[15] H. V. Peiris et ai, "First year Wilkinson Microwave Anisotropy Probe (WMAP) 
observations: Implications for inflation," Astrophys. J. Suppl. 148, 213 (2003), 
astro-ph/0302225. 

[16] J. Preskill, M. B. Wise and F. Wilczek, "Cosmology Of The Invisible Axion," 
Phys. Lett. B 120, 127 (1983); L. F. Abbott and P. Sikivie, "A Cosmological 
Bound On The Invisible Axion," Phys. Lett. B 120, 133 (1983); M. Dine and 
W. Fischler, "The Not-So-Harmless Axion," Phys. Lett. B 120, 137 (1983). 

[17] M. S. Turner, "On The Cosmic And Local Mass Density Of 'Invisible' Axions," 
Phys. Rev. D 33, 889 (1986). 

[18] M. S. Turner, "Windows On The Axion," Phys. Rept. 197, 67 (1990). 



48 



[19] M. Kawasaki, T. Moroi and T. Yanagida, "Can Decaying Particles Raise the 
Upperbound on the Peccei-Quinn Scale?," Phys. Lett. B 383, 313 (1996), hep- 
ph/9510461. 

[20] D. J. Gross, R. D. Pisarski and L. G. Yaffc, "QCD And Instantons At Finite 
Temperature," Rev. Mod. Phys. 53, 43 (1981). 

[21] D. N. Spcrgcl et ai, "First Year Wilkinson Microwave Anisotropy Probe 
(WMAP) Observations: Determination of Cosmological Parameters," astro- 
ph/0302209. 

[22] F. Wilczck, "A model of anthropic reasoning, addressing the dark to ordinary 
matter coincidence," hep-ph/0408167. 

[23] G. Lazarides, R. K. Schaefer, D. Seckel and Q. Shafi, "Dilution Of Cosmological 
Axions By Entropy Production," Nucl. Phys. B 346, 193 (1990). 

[24] M. H. Reno and D. Seckel, "Primordial Nucleosynthesis: The Effects Of Injecting 
Hadrons," Phys. Rev. D 37, 3441 (1988). 

[25] J. E. Kim, "Effects of decay of scalar partner of axion on cosmological bounds of 
axion supermultiplet properties," Phys. Rev. Lett. 67, 3465 (1991); D. H. Lyth, 
"Dilution of cosmological densities by saxino decay," Phys. Rev. D 48, 4523 
(1993), hep-ph/9306293. 

[26] T. Banks, M. Dine and M. Graesser, "Supersymmetry, axions and cosmology," 
hep-ph/0210256. 

[27] M. Dine, W. Fischler and M. Srednicki, "A Simple Solution To The Strong CP 
Problem With A Harmless Axion," Phys. Lett. B 104, 199 (1981). 

[28] K. Rajagopal, M. S. Turner and F. Wilczek, "Cosmological Implications Of Ax- 
inos," Nucl. Phys. B 358, 447 (1991); L. Covi, J. E. Kim and L. Roszkowski, 
"Axinos as cold dark matter," Phys. Rev. Lett. 82, 4180 (1999) hep-ph/9905212. 

[29] T. S. Bunch and P. C. Davies, "Quantum Field Theory In De Sitter Space: 
Renormahzation By Point Splitting," Proc. Roy. Soc. Lond. A 360, 117 (1978). 

[30] L. F. Abbott and M. Wise, "Constraints On Generalized Inflationary Cosmolo- 
gies," Nucl. Phys. B 244, 541 (1984). 

[31] S. M. Leach, A. R. Liddle, J. Martin and D. J. Schwarz, "Cosmological parameter 
estimation and the inflationary cosmology," Phys. Rev. D 66, 023515 (2002), 
astro-ph/0202094. 



49 



[32] A. R. Liddle and D. H. Lyth, "Cosmological inflation and large-scale structure," 
Cambridge, UK: Univ. Pr. (2000) 4OO p 

[33] For a recent review, and references see: M. Kamionkowski, "Inflation at the 
edges," astro-ph/0209273. 

[34] A. R. Cooray, "Polarized CMB: Reionization and primordial tensor modes," 
astro-ph/0311059. 

[35] L. Knox and Y. S. Song, "A limit on the detectability of the energy scale of 
inflation," Phys. Rev. Lett. 89, 011303 (2002), astro-ph/0202286. 

[36] E. W. Kolb and M. S. Turner, "The Early Universe," Redwood City, USA: 
Addison- Wesley (1990) 5^7 p. (Frontiers in physics, 69). 

[37] R. K. Sachs and A. M. Wolfe, "Perturbations Of A Cosmological Model And An- 
gular Variations Of The Microwave Background," Astrophys. J. 147, 73 (1967). 

[38] J. R. Bond and G. Efstathiou, "The Statistics Of Cosmic Background Radiation 
Fluctuations," Mon. Not. Roy. Astron. Soc. 226, 655 (1987). 

[39] U. Seljak and M. Zaldarriaga, http://www.cmbfast.org, 
http: / /lambda.gsfc.nasa.gov/toolbox/tb_cmbfast_ov.cfm 

[40] C. L. Bennett et al, "4- Year COBE DMR Cosmic Microwave Background Ob- 
servations: Maps and Basic Results," Astrophys. J. 464, LI (1996), astro- 
ph/9601067. 

[41] P. Crotty, J. Garcia-BeUido, J. Lesgourgues and A. Riazuelo, "Bounds on isocur- 
vature perturbations from CMB and LSS data," astro-ph/0306286. 

[42] E. Komatsu et al, "First Year Wilkinson Microwave Anisotropy Probe (WMAP) 
Observations: Tests of Gaussianity," astro-ph/0302223. 

[43] M. Dine, L. Randall and S. Thomas, "Supersymmetry breaking in the early 
universe," Phys. Rev. Lett. 75, 398 (1995), hep-ph/9503303. 

[44] R. A. Battye and E. P. S. Shellard, "Axion string constraints," Phys. Rev. Lett. 
73, 2954 (1994) [Erratum-ibid. 76, 2203 (1996)], astro-ph/9403018. 

[45] G. G. Raffelt, "Astrophysical Methods To Constrain Axions And Other Novel 
Particle Phenomena," Phys. Rept. 198 (1990) 1. 

[46] M. Dine and A. Anisimov, "Is there a Peccei-Quinn phase transition?," hep- 
ph/0405256. 



50 



[47] S. Dimopoulos and L. J. Hall, "Inflation And Invisible Axions," Phys. Rev. Lett. 
60, 1899 (1988). 

[48] E. Witten, "Dimensional Reduction Of Superstring Models," Phys. Lett. B 155, 
151 (1985). 

[49] K. Choi and J. E. Kim, "Harmful Axions In Superstring Models," Phys. Lett. 
B 154, 393 (1985) [Erratum-ibid. 156B, 452 (1985)]; K. Choi and J. E. Kim, 
"Compactification And Axions In E(8) X E(8)-Prime Superstring Models," Phys. 
Lett. B 165, 71 (1985). 

[50] M. Srcdnicki, "Axion Couplings To Matter. 1. CP Conserving Parts," Nucl. Phys. 
B 260, 689 (1985). 



51 



